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1 Introduction 



o 
o 

(N 

$H ' Abstract. A perturbative approach to quantum field theory involves 

C^' the computation of loop integrals, as soon as one goes beyond the leading 

■^C , term in the perturbative expansion. First I review standard techniques 

for the computation of loop integrals. In a second part I discuss more 

advanced algorithms. For these algorithms algebraic methods play an 

important role. A special section is devoted to multiple polylogarithms. 

►^ ' I tried to make these notes self-contained and accessible both to 

QQ ' physicists and mathematicians. 

o 
^ 

o 

^D I In this lecture I will discuss techniques for the computation of loop integrals, 

which occur in perturbative calculations in quantum field theory. But before em- 
C^hI barking onto a journey of integration and special function theory, it is worth recall- 

ing the motivation for such an effort. 
^ I High-energy physics is successfully described by the Standard Model. The term 

r^ • "Standard Model" has become a synonym for a quantum field theory based on the 

gauge group SU{3) (E) SU{2) (g) C/(l). At high energies all coupHng constants are 
. small and perturbation theory is a valuable tool to obtain predictions from the 

rS ' theory. For the Standard Model there are three coupling constants, gi, (72 and 173, 

S i corresponding to the gauge groups C/(l), SU{2) and SU{3), respectively. As all 

methods which will be discussed below, do not depend on the specific nature of 
these gauge groups and are even applicable to extensions of the Standard Model 
(like supersymmetry) , I will just talk about a single expansion in a single coupling 
constant. All observable quantities are taken as a power series expansion in the 
coupling constant, and calculated order by order in perturbation theory. 

Over the years particle physics has become a field where precision measurements 
have become possible. Of course, the increase in experimental precision has to be 
matched with more accurate calculations from the theoretical side. This is the 
"raison d'etre" for loop calculations: A higher accuracy is reached by including 
more terms in the perturbative expansion. There is even an additional "bonus" we 
get from loop calculations: Inside the loops we have to take into account all particles 
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which could possibly circle there, even the ones which are too heavy to be produced 
directly in an experiment. Therefore loop calculations in combination with precision 
measurements allow us to extend the range of sensitivity of experiments from the 
region which is directly accessible towards the range of heavier particles which 
manifest themselves only through quantum corrections. As an example, the mass 
of top quark has been predicted before the discovery of the top quark from the loop 
corrections to electro-weak precision experiments. The same experiments predict 
currently a range for the mass of the yet undiscovered Higgs boson. 

It is generally believed that a perturbative series is only an asymptotic series, 
which will diverge, if more and more terms beyond a certain order are included. 
However this shall be of no concern to us here. We content ourselves to the first few 
terms in the perturbative expansion with the implicit assumption, that the point 
where the power series starts to diverge is far beyond our computational abilities. In 
fact, our computational abilities are rather limited. The complexity of a calculation 
increases obviously with the number of loops, but also with the number of external 
particles or the number of non-zero internal masses associated to propagators. To 
give an idea of the state of the art, specific quantities which are just pure numbers 
have been computed up to an impressive fourth or third order. Examples are the 
calculation of the 4-loop contribution to the QCD /?- function [1], the calculation 
of the anomalous magnetic moment of the electron up to three loops [2] , and the 
calculation of the ratio 

cr(e+e^ -^ hadrons) 

R = — _ — — — (1.1) 

cr(e+e -^ ^i+fi ) 

of the total cross section for hadron production to the total cross section for the pro- 
duction of a n^fi~ pair in electron-positron annihilation to order O [g^j (also involv- 
ing a three loop calculation) [3]. Quantities which depend on a single variable are 
known at the best to the third order. Outstanding examples are the computation of 
the three-loop Altarelli-Parisi splitting functions [4,5] or the calculation of the two- 
loop amplitudes for the most interesting 2—^2 processes [6-16]. The complexity of 
a two-loop computation increases, if the result depends on more than one variable. 
An example for a two-loop calculation whose result depends on two variables is the 
computation of the two-loop amplitudes for e~^e~ —^ 3 jets [17-19]. But in general, 
if more than one variable is involved, we have to content ourselves with next-to- 
leading order calculations. An example for the state of the art is here the compu- 
tation of the electro- weak corrections to the process e+e^ -^ 4 fermions [20,21]. 

From a mathematical point of view it is an interesting question to ask which 
type of functions or constants appear in the result of a particular calculation. For 
one-loop calculations it is known that the result can be expressed in terms of log- 
arithms and dilogarithms. Transcendental constants which do occur are just pro- 
portional to TT^, or phrased differently, proportional to C2 = 7r^/6. Going to higher 
loops one finds the following: In results which arc just numbers higher zeta values 
do occur. In results which depend on one variable harmonic polylogarithms arise 
and in results which depend on more than one variable multiple polylogarithms 
occur. The multiple polylogarithms have nice algebraic properties and contain as 
special cases (multiple) zeta values and the harmonic polylogarithms. They will be 
discussed in detail in this lecture. 

This paper is organised as follows: In the next section I will discuss basic tech- 
niques, which allow us to exchange the integrals over the loop momenta against 
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integrals over Feynman or Schwinger parameters. Sect. |2| shows how tensor inte- 
grals, e.g. integrals where the loop momentum occurs also in the numerator, can be 
reduced to scalar integrals. In sect. 0)1 discuss how the Feynman parametrisation for 
a generic scalar Z-loop integral can be read off directly from the underlying Feynman 
graph. The simplest, but most important loop integrals are the one-loop integrals, 
whose theory is presented in sect.O The main results are that all one- loop integrals 
can be reduced to one-loop integrals with no more than four external legs. Fur- 
thermore, when computed within dimensional regularisation up to the finite part, 
the only occurring transcendental functions are the logarithm and the dilogarithm. 
In sect.inil continue with the general case of an /-loop integral and present several 
advanced methods for the calculation of the Feynman parameter integrals. As al- 
ready mentioned, the multiple polylogarithms are the class of functions, to which 
the logarithm and the dilogarithm naturally extend. They are discussed in detail 
in sect |3 Unfortunately we do not yet have a complete theory for loop integrals in 
quantum theory. There are many unknown integrals and a variety of conjectures 
related to them. Sect.|Hltries to give an outlook towards open questions and future 
directions. Finally, sect. |21 contains a summary. 

2 Basic techniques 

To set the scene, let us start with a brief summary how observables are cal- 
culated in perturbative quantum field theory. Let us assume that we are inter- 
ested in a process with two incoming elementary particles (like an electron and 
a positron) and n outgoing particles. We study an observable O {pi,p2;ki, ...,kn) 
which depends on the four-momenta pi and p2 of the two incoming particles and 
the four-momenta kj of the outgoing particles and we are interested in the value of 
O integrated over all possible final state momenta kj. O {pi,p2;ki, ..., fc„) depends 
on the experimental set-up and can be an arbitrary complicated function of the 
four-momenta. In the simplest case this function is just a constant equal to one, 
corresponding to the situation where we count every event with n particles in the 
final state. In more realistic situations one takes for example into account that it is 
not possible to detect particles close to the beam pipe. The function O would then 
be zero in these regions of phase space. 

The expectation value for the observable O is given by 

(2.1) 

where l/8/(pi +P2)^ is a normalisation factor taking into account the incoming flux 
and averages over the spins of the incoming particles. The phase space measure is 
given by 

1 " J37, . / " \ 

There is a symmetry factor 1/Njl for each type of particle, which occurs more than 
once in the final state. Ei is the energy of particle i: 



E, = \^kf+ml (2.3) 
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As the integrand can be a rather compHcated function, the phase space integral is 
usuaUy performed numerically by Monte Carlo integration. 

The most important ingredient of formula 1)2. l|l is the matrix element An for 
the process under consideration. This quantity is also called the amplitude or the 
Green function of the process. The amplitude has the perturbative expansion 

An = 9-{A^^^+9'Ail^+gU^:^+g'A^^^ + ...) (2.4) 

To the coefficient Ah, contribute Feynman graphs with I loops and (n + 2) external 
legs. The recipe for the computation of An is as follows: Draw first all Feynman 
diagrams with the given number of external particles and I loops. Then translate 
each graph into a mathematical formula with the help of the Feynman rules. An 
is then given as the sum of all these terms. 



2.1 Feynman rules. Feynman rules allow us to translate a Feynman graph 
into a mathematical formula. These rules are derived from the fundamental La- 
grange density of the theory, but for our purposes it is sufficient to accept them as 
a starting point. The most important ingredients are internal propagators, vertices 
and external lines. For example, the rules for the propagators of a fermion or a 
gauge boson read 

„ . . ^+m 

i'ermion: ' = i- 



p^ — mP- + i5 ' 
Gauge boson: ^WW = _1_ (^g^, _ (i _ ^^^ . (2.5) 

Here p and k are the momenta of the fermion and the boson, respectively, m 
is the mass of the fermion. ^ = Pfil^ is a short-hand notation for the con- 
traction of the momentum with the Dirac matrices. The metric tensor is de- 
noted by g^i, and the convention adopted here is to take the metric tensor as 
9ij.u — diag(l, — 1, — 1, — 1). The propagator would have a pole for p^ = m? , or 
phrased differently E — ± yjfp + m? . When integrating over E, the integration 
contour has to be deformed to avoid these two poles. Causality dictates into which 
directions the contour has to be deformed. The pole on the negative real axis is 
avoided by escaping into the lower complex half-plane, the pole at the positive real 
axis is avoided by a deformation into the upper complex half-plane. Feynman in- 
vented the trick to add a small imaginary part i5 to the denominator, which keeps 
track of the directions into which the contour has to be deformed. I will usually 
suppress the Z(5-term in order to keep the notation compact. 

The gauge boson propagator depends on an additional variable ^, called the 
gauge-fixing parameter. This parameter can be chosen arbitrarily, but one must 
use the same choice for all diagrams contributing to an gauge-invariant subset of 
the amplitude. It also shows that there is no point in speaking about the "value" 
of an individual graph, as this value depends in general on our choice of the gauge 
parameter. Only in the sum of graphs over a gauge-invariant set does this depen- 
dency on the gauge parameter cancel. The most economical choice is ^ = 1, which 
is called Feynman gauge. 
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P4 / ^ P3 

Figure 1 A one-loop Feynman diagram contributing to tiie process e'^e~ — > qgq. 

As a typical example for an interaction vertex let us look at the vertex involving 
a fermion pair and a gauge boson: 



igj^. (2.6) 



Here, g is the coupling constant and 7^ denotes the Dirac matrices. At each vertex, 
we have momentum conservation: The sum of the incoming momenta equals the 
sum of the outgoing momenta. 

To each external line we have to associate a factor, which describes the polar- 
isation of the corresponding particle: There is a polarisation vector e''(fc) for each 
external gauge boson and a spinor m(p), u(p), v{p) or v{p) for each external fermion. 

Furthermore there are a few additional rules: First of all, there is an integration 

'"' (2.7) 




(27r)4 

for each loop. Secondly, each closed fermion loop gets an extra factor of (—1). 
Finally, each diagram gets multiplied by a symmetry factor 1/ S, where S is the 
order of the permutation group of the internal lines and vertices leaving the diagram 
unchanged when the external lines are fixed. 

Having stated the Feynman rules, let us look at an example: Fig. ^ shows a 
Feynman diagram contributing to the one-loop corrections for the process e~^e~ — > 
qgq. At high energies we can ignore the masses of the electron and the light quarks. 
From the Feynman rules one obtains for this diagram (ignoring coupling and colour 
prefactors): 

~v{pa)-i^u{pz)^ j j^-^u{pi)^{p2)^lu^lt.-^Yv{pi). (2.8) 

Here, pi2 = Pi + P2, P123 ^ Pi + P2 + P3, h = fci - P12, fcs == ^2 - Ps- Further 
^(^2) = 7re^(P2), where £^(^2) is the polarisation vector of the outgoing gluon. 
All external momenta are assumed to be massless: pf = for i = 1..5. We can 
reorganise this formula into a part, which depends on the loop integration and a 
part, which docs not. The loop integral to be calculated reads: 

a Hi n^n^ 
while the remainder, which is independent of the loop integration is given by 

- V{p4:)j''u{p5)— —U{pi)^{p2)^12l,.lpltj.ja7''v{p3). (2.10) 

P123P12 
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The loop integral in eq. (|2.9() contains in the denominator three propagator factors 
and in the numerator two factors of the loop momentum. We call a loop integral, 
in which the loop momentum occurs also in the numerator a "tensor integral" . 
A loop integral, in which the numerator is independent of the loop momentum is 
called a "scalar integral" . The basic strategy, which will be discussed in detail in 
sect. 121 consists in reducing tensor integrals to scalar integrals. The scalar integral 
associated to cq. (12.911 reads 



d^ki 



l^TTI rC-i Kn rvr) 



noT^- (2-11) 



2.2 Regularisation. Before we start with the actual calculation of loop inte- 
grals, I should mention one complication: Loop integrals are often divergent ! Let 
us first look at the simple example of a scalar two-point one-loop integral with zero 
external momentum: 



p = 




d^k 1 
(27r)4 (fc2)2 

OO OO 



(47r)2 J k^ (47r)2 J x 



This integral diverges at fc^ ^ oo as well as at A:^ — > 0. The former divergence is 
called ultraviolet divergence, the later is called infrared divergence. Any quantity, 
which is given by a divergent integral, is of course an ill-defined quantity. There- 
fore the first step is to make these integrals well-defined by introducing a regulator. 
There are several possibilities how this can be done, but the method of dimensional 
regularisation [22-24] has almost become a standard, as the calculations in this 
regularisation scheme turn out to be the simplest. Within dimensional regulari- 
sation one replaces the four-dimensional integral over the loop momentum by an 
£)-dimensional integral, where D is now an additional parameter, which can be a 
non-integer or even a complex number. We consider the result of the integration 
as a function of D and we are interested in the behaviour of this function as D 
approaches 4. The D-dimensional integration still fulfils the standard laws for in- 
tegration, like linearity, translation invariance and scaling behaviour [25,26]. If / 
and g are two functions, and if a and h are two constants, linearity states that 



d^k {af{k) + bg{k)) = a d^kf{k) + b / d"kg{k). (2.13) 

Translation invariance requires that 

fd"kf{k+p) = fd"kf{k). (2.14) 

for any vector p. 

The scaling law states that 

d^kf{\k) = A"-" f d^kf{k). (2.15) 
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The D-dimensional integral has also a rotation invariance: 

d^kJ{Kk) = f d^kf{k), (2.16) 

where A is an element of the Lorentz group 5*0(1, 1? — 1) of the I?-diniensional 
vector-space. Here we assumed that the D-dimensional vector-space has the metric 
diag(-|-l, — 1, — 1, — 1, ...). The integral measure is normalised such that it agrees 
with the result for the integration of a Gaussian function for all integer values D: 

D_ 

i^A:exp(afc2) = ^ (-) ' • (2-17) 

We will further assume that we can always decompose any vector into a 4-dimcnsional 
part and a, [D — 4)-dimensional part 

^(D) ^ ^(4) + ^(13-4)' (2-18) 

and that the 4-dimensional and [D — 4)-dimcnsional subspaces are orthogonal to 
each other: 

%) • fc(D-4) - 0. (2.19) 

If D is an integer greater than 4, this is obvious. We postulate that these relations 
arc true for any value of D. One can think of the underlying vector-space as a space 
of infinite dimension, where the integral measure mimics the one in D dimensions. 

In practice we will always arrange things such that every function we integrate 
over D dimensions is rotational invariant, e.g. is a function of k^ . In this case 
the integration over the [D — 1) angles is trivial and can be expressed in a closed 
form as a function of D. Let us assume that we have an integral, which has a 
UV-divergence, but no IR-divergences. Let us further assume that this integral 
would diverge logarithmically, if we would use a cut-off rcgularisation instead of 
dimensional rcgularisation. It turns out that this integral will be convergent if the 
real part of D is smaller than 4. Therefore we may compute this integral under 
the assumption that Re(£') < 4 and we will obtain as a result a function of D. 
This function can be analytically continued to the whole complex plane. We are 
mainly interested in what happens close to the point D = A. For an ultraviolet 
divergent one-loop integral we will find that the analytically continued result will 
exhibit a pole at D = 4. It should be mentioned that there are also integrals which 
are quadratically divergent, if a cut-off regulator is used. In this case we can repeat 
the argumentation above with the replacement Re(£') < 2. 

Similarly, we can consider an IR-divergent integral, which has no UV-divergence. 
This integral will be convergent if Re(-D) > 4. Again, wc can compute the inte- 
gral in this domain and continue the result to D = 4. Here we find that each 
IR-divergent loop integral can lead to a double pole a.i D — A. 

We will use dimensional rcgularisation to regulate both the ultraviolet and 
infrared divergences. The attentative reader may ask how this goes together, as 
we argued above that UV-divergences require Re(£') < 4 or even Re(£') < 2, 
whereas IR-divergences are regulated by Re(Z?) > 4. Suppose for the moment that 
we use dimensional rcgularisation just for the UV-divergences and that we use a 
second regulator for the IR-divergences. For the IR-divergences we could keep all 
external momenta off-shell, or introduce small masses for all massless particles or 
even raise the original propagators to some power v. The exact implementation 
of this regulator is not important, as long as the IR-divergences are screened by 
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this procedure. We then perform the loop integration in the domain where the 
integral is UV-convergent. We obtain a result, which we can analytically continue 
to the whole complex Z3-plane, in particular to Re(£') > 4. There we can remove 
the additional regulator and the IR-divergences are now regulated by dimensional 
regularisation. Then the infrared divergences will also show up as poles at D — A. 

In summary, within dimensional regularisation the initial divergences show up 
as poles in the complex Z3-plane at the point Z) = 4. What shall we do with these 
poles ? The answer has to come from physics and we distinguish again the case 
of UV-divergences and IR-divergences. The UV-divergences are removed through 
renormalisation. On the level of Feynman diagrams we can associate to any UV- 
divergent part a counter-term, which has exactly the same pole term at _D = 4, but 
with an opposite sign, such that in the sum the two pole terms cancel. 

As far as infrared-divergences are concerned we first note that any detector 
has a finite resolution. Therefore two particles which are sufficiently close to each 
other in phase space will be detected as one particle. Now let us look again at 
eqs. (|2.1|) and (|2.4|) . The next-to- leading order term will receive contributions 
from the interference term of the one-loop amplitude Ah, with the leading-order 
amplitude An , both with n final state particles. This contribution is of order 
g'2n+2^ Of the same order is the square of the leading-order amplitude A^^^i with 
(n+l) final state particles. This contribution we have to take into account whenever 
our detector resolves only n particles. It turns out that the phase space integration 
over the regions where one or more particles become unresolved is also divergent, 
and, when performed in D dimensions, leads to poles with the opposite sign as the 
one encountered in the loop amplitudes. Therefore the sum of the two contributions 
is finite. The Kinoshita-Lee-Nauenberg theorem [27,28] guarantees that all infrared 
divergences cancel, when summed over all degenerate physical states. To make this 
cancellation happen in practice requires usually quite some work, as the different 
contributions live on phase spaces of different dimensions. 

To summarise we are interested into loop integrals regulated by dimensional 
regularisation. As a result we seek the Laurent expansion around D = A. It is 
common practice to parameterise the deviation of D from 4 by 

D = 4-2e. (2.20) 

Divergent loop integrals will therefore have poles in 1/e. In an Z-loop integral ultra- 
violet divergences will lead to poles 1/e' at the worst, whereas infrared divergences 
can lead to poles up to 1/e^'. We will also encounter integrals, where the dimension 
is shifted by units of two. In these cases we often write 

D = 2m~2e, (2.21) 

where m is an integer, and we are again interested in the Laurent series in e. 

2.3 Dimensional regularisation schemes. Let us look again at the exam- 
ple in eqs. H2.8|l to H2.1Q() . We separated the loop integral from the remainder in 
eq. (12.10(1 . which is independent of the loop integration. In this remainder the 
following string of Dirac matrices occurs: 

l.lpl^lal". (2.22) 
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If we anti-commute the first Dirac matrix, we can achieve that the two Dirac ma- 
trices with index v are next to each other: 

7.7"- (2.23) 

In four dimensions this equals 4 times the unit matrix. What is the value within 
dimensional regularisation ? The answer depends on how we treat the Dirac algebra. 
Does the Dirac algebra remain in four dimensions or do we also continue the Dirac 
algebra to D dimensions ? There are several schemes on the market which treat 
this issue differently. To discuss these schemes it is best to look how they treat 
the momenta and the polarisation vectors of observed and unobserved particles. 
Unobserved particles are particles circulating inside loops or emitted particles not 
resolved within a given detector resolution. The most commonly used schemes are 
the conventional dimensional regularisation scheme (CDR) [26], where all momenta 
and all polarisation vectors are taken to be in D dimensions, the 't Hooft-Veltman 
scheme (HV) [22,29], where the momenta and the helicities of the unobserved 
particles are D dimensional, whereas the momenta and the helicities of the observed 
particles are 4 dimensional, and the four-dimensional helicity scheme (FD) [30-32], 
where all polarisation vectors are kept in four dimensions, as well as the momenta of 
the observed particles. Only the momenta of the unobserved particles are continued 
to D dimensions. 

The conventional scheme is mostly used for an analytical calculation of the 
interference of a one-loop amplitude with the Born amplitude by using polarisa- 
tion sums corresponding to D dimensions. For the calculation of one-loop helicity 
amplitudes the 't Hooft-Veltman scheme and the four-dimensional helicity scheme 
are possible choices. All schemes have in common, that the propagators appearing 
in the denominator of the loop-integrals are continued to D dimensions. They dif- 
fer how they treat the algebraic part in the numerator. In the 't Hooft-Veltman 
scheme the algebraic part is treated in D dimensions, whereas in the FD scheme 
the algebraic part is treated in four dimensions. It is possible to relate results 
obtained in one scheme to another scheme, using simple and universal transition 
formulae [33-35]. 

Therefore, if we return to the example above, we have 

I, _ ( D ■ 1, in the CDR and HV scheme, , , 

^•"^ ^ \ 4 • 1, in the FD scheme. ^ ' 

2.4 Loop integration in D dimensions. In this section I will discuss how 
to perform the D-dimensional loop integrals. It would be more correct to say that 
we exchange them for some parameter integrals. Our starting point is a one-loop 
integral with n external legs: 

^"v'7^ (_p,(-p.)..^(-p,.) - ''■''> 

where the propagators are of the form 

P, = ik-J2pA ~^1 (2-26) 

and pj are the external momenta. The small imaginary parts i6 are not written 
explicitly. In eq. (|2.25|l there are some overall factors, which I inserted for con- 
venience: /i is an arbitrary mass scale and the factor /i^"^ ensures that the mass 
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dimension of eq. H2.25() is an integer. The factor e'^^'^ avoids a proliferation of 
Euler's constant ^e in the final result. The integral measure is now d^k/{i7r^'^) 
instead of dPk/{2Ti)^ , and each propagator is multiplied by (—1). The reason for 
doing this is that the final result will be simpler. 

In order to perform the momentum integration we proceed by the following 
steps: 

1. Feynman or Schwinger parametrisation. 

2. Shift of the loop momentum to complete the square, such that the integrand 
depends only on fc^. 

3. Wick rotation. 

4. Introduction of generalised spherical coordinates. 

5. The angular integration is trivial. Using the definitions of Euler's Gamma 
and Beta functions, the radial integration can be performed. 

6. This leaves only the non-trivial integration over the Feynman parameters. 

Although I discuss here only one-loop integrals, the methods presented in this 
section are rather general and can be applied iteratively to /-loop integrals. 

2.4.1 Feynman and Schwinger parameterisation. As already discussed above, 
the only functions we really want to integrate over D dimensions are the ones which 
depend on the loop momentum only through k^. The integrand in eq. (|2.25|) is not 
yet in such a form. To bring the integrand into this form, we first convert the 
product of propagators into a sum. To do this, there are two techniques, one due 
to Feynman, the other one due to Schwinger. Let me start with the Feynman 
parameter technique. In its full generality it is also applicable to cases, where each 
factor in the denominator is raised to some power v. The formula reads: 






1=1 



i=l 



(2.27) 



The proof of this formula can be found in many text books and is not repeated 
here. The price we have to pay for converting the product into a sum are (n — 1) 
additional integrations. Let us look at the example from eq. (|2.11|) : 

' 2ld./fd./fd., /^'-^^~T~^t. (2.28) 



(-fc?)(-fci)Mi) J J J (-x,ki-x2ki~x,kif 

P ^ ^ z -^ 6/ 

1 l-Kl 

= 2 / dxi / dx2 







{ — Xik\ — X2fc| — (1 — Xi — X2)k'l) 



An alternative to Feynman parameters are Schwinger parameters. Here each prop- 
agator is rewritten as 



1 







CSO 

l^/.„'-exp,.P,. ,2.2,) 
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Im kn 




Re kn 



Figure 2 Integration contour for the Wick rotation. The little circles along 
the real axis exclude the poles. 



Therefore we obtain for our example 
1 



{-kl)i-kl)i-kl) 



/ dxi / dx2 / dx^ exp (xifcj + a;2A;| + xafcg) . (2.30) 



2.4.2 Shift of the integration variable. We can now complete the square and 
shift the loop momentum, such that the integrand becomes a function of k^ . This 
is best discussed by an example. We consider again eq. (|2.28l) . With fc2 = fci — pi2 
and fcs = fc2 ~ P3 we have 

- xikl - X2kl - x^kl = - {ki - X2P12 - X3P123) -a;iX2Si2 -a;ia;3Si23,(2.31) 

where S12 = (pi +^2)^ and S123 = {pi + P2 +^3)^- We can now define 

k[ = fci - a;2Pi2 - 2:3^123 (2.32) 

and using translational invariancc our loop integral becomes 

d^ki 1 



111 
^ ' ■ D/\ / ^^1 / ^^2 / dX3 



6{1 -Xx- X2 ~ X-i) 



{~K 



xiX2S\2 ~ a;ia;3Si23 



(2.33) 



The integrand is now a function of k'^ ■ 

2.4.3 Wick rotation. Having succeeded to rewrite the integrand as a function 
of fc^, we then perform a Wick rotation, which transforms Minkowski space into an 
Euclidean space. Remember, that k'^ written out in components in D- dimensional 
Minkowski space reads 



ftp 



kl 



^2 



fcg 



(2.34) 



(Here kj denotes the j'-th component of the vector fc, in contrast to the previous 
section, where we used the subscript to label different vectors kj. It should be 
clear from the context what is meant.) Furthermore, when integrating over fcg, we 
encounter poles which are avoided by Feynman's i(5-prescription. In the complex 
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A:o-plane we consider the integration contour shown in fig. [21 Since the contour does 
not enclose any poles, the integral along the complete contour is zero: 

dkofiko) = 0. (2.35) 

If the quarter-circles at infinity give a vanishing contribution (it can be shown that 
this is the case) we obtain 

— zoo 

dkofih) = ~ I dkofiko). (2.36) 

We now make the following change of variables: 

fco = iKo, 

kj = Kj, forl<j<D-l. (2.37) 



As a consequence we have 



fc2 = -i^2^ 



d'^k = id^K, (2.38) 

where K^ is now given with Euclidean signature: 

K"^ = Kl + Kl+Kl + Kl + ... (2.39) 

Combining eq. H2.36|l with eq. I|2.37|l we obtain for the integration of a function 
/(fc^) over D dimensions 

rl^h f d^ K 

-J{-k^) - / -^f{K% (2.40) 



iT,D/2''^ > J ^D/2- 

whenever there are no poles inside the contour of fig. [5] and the arcs at infinity 
give a vanishing contribution. The integral on the r.h.s. is now over Z?-dimensional 
Euclidean space. Eq. (|2.4U|I justifies our conventions, to introduce a factor i in the 
denominator and a minus sign for each propagator in eq. H2.25|l . These conventions 
are just such that after Wick rotation wc have simple formulae. 

2.4.4 Generalised spherical coordinates. We now have an integral over D-dimen- 
sional Euclidean space, where the integrand depends only on K^. It is therefore 
natural to introduce spherical coordinates. In D dimensions they are given by 

Ko — K cos 9i , 

Ki = iiT sin 01 cos 6*2 , 

Kd-2 = i^sin0i...sin6'_D_2COS0£)_i, 

Kd^i = Ks\n6i...sin9D-2S\n9D~i- (2.41) 

In D dimensions we have one radial variable K , D — 2 polar angles 9j (with 1 < 
j < D — 2) and one azimuthal angle 0d-i. The measure becomes 

d^K = K^-^dKdflD, (2.42) 

where 

D-l 

d^D = W&m'^'^^'e.de,. (2.43) 
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Integration over the angles yields 

f d^D - f dOisin^'-^ei... f deD-2sm9D-2 f d9D-i = ^y^, (2.44) 



where T{x) is Euler's Gamma fmiction. Note that the integration on the l.h.s of 
eq. (|2.44ll is defined for any natural number D, whereas the result on the r.h.s is 
an analytic function of D, which can be continued to any complex value. 

2.4.5 Euler's Gamma and Beta function. It is now the appropriate place to 
introduce two special functions, Euler's Gamma and Beta function, which are used 
within dimensional regularisation to continue the results from integer D towards 
non-integer values. The Gamma function is defined for Re(a:) > by 



(2.45) 



r(a;) = / e-H'^-^dt. 
Jo 

It fulfils the functional equation 

r(a; + l) = xT{x). (2.46) 

For positive integers n it takes the values 

r(n + l) = n! = l-2-3-...-n. (2.47) 

At X = 1/2 it has the value 

rfi) = V^, (2-48) 



2^ 
which can also be inferred from the relation 



r(a;)r(l-2;) = . (2.49) 

SmTTX 

For integers n we have the reflection identity 

Y{x-n) r(i-x) 

^X^ - (-1) T{l-x + ny (2-^^) 

The Gamma function Tix) has poles located on the negative real axis at x — 
0, —1, —2, .... Quite often we will need the expansion around these poles. This can 
be obtained from the expansion around x = 1 and the functional equation. The 
expansion around e = 1 reads 

r(l+e) = exp I -7££ + ^ ^-^Cn£" 1 , (2.51) 

where 'ys is Euler's constant 



7b 
and Cn is given by 



= lim I V- -Inn =0.5772156649... (2.52) 



oo ^ 

E-^- (2-53) 



For example we obtain for the Laurent expansion around e = 



r(£) - i - 7B + 0(e). (2.54) 

e 
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Euler's Beta function is defined for Re(a;) > and Re(y) > by 

1 
B{x,y) = ft''-\l-t)y-^dt, (2.55) 



or equivalently by 

oo 

B(x,y) ^ I -^———dt. (2.56) 


The Beta function can be expressed in terms of Gamma functions: 

''<-) ^ Wfy <'■"> 

2.4.6 Result for the momentum integration. We are now in a position to per- 
form the integration over the loop momentum. Let us discuss again the example 
from eq. (|2.33|) . After Wick rotation we have 



(2.58) 

Introducing spherical coordinates and performing the angular integration this be- 
comes 

CO -P-2 

2 f ^ f ^ 5(1 - xi - X2 - x:i) (K^) " 
r - ^ / dK^ / d^x- --^ ^^-. (2.59) 



^ (y) J J {K'^ ~ X1X2S12 ~ xiX'iSi 



23 J 



For the radial integration we have after the substitution t = K'^/{—xiX2Si2 
a;ia;3Si23) 



D-2 



D-2 



dK'^- — -3- = (-xia;2Si2 - a;ia;3Si23) ^ ^ l dt 



[K^ - X1X2S12- X1X3S123) J (1 + i) 

u 

(2.60) 
The remaining integral is just the second definition of Euler's Beta function 

'\ t^ r(f)r(3-f) , , 

dt =■ = ^^^ \ 2J.. (2.61 

{1 + tf r(3) ^ ' 



Putting everything together and setting D = 4 — 2e we obtain 
d^ki 1 

r (1 + e) / d^x (5(1 - xi - X2 - ^3) x^^^" {-X2S12 - X3Si23y^~^ 



(2.62) 



Therefore we succeeded in performing the integration over the loop momentum k 
at the expense of introducing a two-fold integral over the Feynman parameters. 

As the steps discussed above always occur in any loop integration we can com- 
bine them into a master formula. IfU and J- are functions, which are independent 



The Art of Computing Loop Integrals 15 

of the loop niomentum, we have for the integration over Minkowski space with 
dimension D ~ 2m — 2e: 



I 






The functions U and J-^ depend usually on the Feynnian parameters and the external 
momenta and are obtained after Feynman parametrisation from completing the 
square. In eq. H2.63|l we allowed additional powers (— A:^)° of the loop momentum 
in the numerator. This is a slight generalisation and will be useful later. Here we 
observe that the dependency of the result on a, apart from a factor T(m + a — 
e)/T{m — e), occurs only in the combination m + a — e = D/2-\-a. Therefore adding 
a power of (— fc^) to the numerator is almost equivalent to consider the integral 
without this power in dimensions D + 2. 

There is one more generalisation: Sometimes it is convenient to decompose k^ 
into a (27TT,)-dimensional piece and a remainder: 

^{D) = ^(2m) + ^(-2e)- (2-64) 

If D is an integer greater than 2m we have 

h^ — i-^ _ jl2 _ _ 1,2 

'^(2m) ^ '^0 '*'! ••• '^2m-n 

fc(_2e) = -4„---4-i- (2.65) 

We also need loop integrals where additional powers of {—kjr^^-.) appear in the 
numerator. These are related to integrals in higher dimensions as follows: 

^2m-2ej^ Y(r — e) r (^2m+2i — 2e^ 

iT^m-e (^^(-2£))'/(fc('2m)'^(-2e)) = ris) J i7r™+^-e ^ ^^2va^ ^ -2e) ■ 

(2.66) 

Here, f{k^2m)'^f-2e)) ^^ ^ function which depends on ^2™, fc2m+i7 •■■, ^d-i only 
through k?_^^y The dependency on fco, fci, ..., A:2m-i is not constrained. 
Finally it is worth noting that 

d^'^-^'k , „.a _ r (-l)"r(a+l), ifm + a-£ = 



^ \ 0, otherwise. ' 

The non-zero value for m + a — £ — can be verified by expanding eq. H2.17JI into 
a power series in a. 

3 Tensor integrals 

In the previous section we considered scalar integrals and integrals where the 
functional dependence on the loop momentum of the numerator of the integrand 
is particular simple, like for example through factors k^ or fc? 2 v ^^ ^^ recall the 
example discussed in eqs. (|2.8|) and (|2.9|l we see that more general tensor structures 
occur. In this section I discuss systematic algorithms for the reduction of tensor 
integrals to scalar integrals. In sect. 13. il l discuss the Passarino-Veltman reduction 
technique, which historically was the first systematic procedure and can be applied 
to one-loop integrals. In sect. 13.21 1 present a more general reduction method, which 
applies to arbitrary Z-loop integrals. 
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3.1 Passarino - Veltman reduction. For one-loop integrals a systematic 
algorithm has been first worked out by Passarino and Veltman [36]. Let us first 
introduce a notation for one-loop tensor integrals. For integrals with one, two or 
three external legs we write 






s^E 2e I "• "' 1, fe^, /C^/C,y 



ITT 



^/2 (-fc2 + mj){-{k - pY + ml) ' 



Co,^,^i/(pi,P2,mi,m2,m3) 



p<^iE ij^s I "' ^ 1, fc^, k^k y 



2\' 



j7r-D/2 (^_}.2 _|_ ml){-{k - pi)2 -f ml){-{k - pi - P2Y + ml) 

with an obvious generalisation towards more external legs and higher rank tensor 
integrals. The notation implies that the numerator for Bq equals 1, for B^ the nu- 
merator equals k^ and for B^j^^, the numerator equals k^ki,. The reduction technique 
according to Passarino and Veltman uses the fact that due to Lorentz symmetry the 
result can only depend on tensor structures which can be build from the external 
momenta p^ and the metric tensor g^^'^ . We therefore write the tensor integrals in 
the most general form in terms of form factors times external momenta and/or the 
metric tensor. For example 

B^" = p^Bi, 

C^-- = pWC21+P^P2'C22 + KP2+KP2)C23+5^''C24. (3.2) 

One then solves for the form factors Bi, 521, ^22, C'li, etc. by first contracting 
both sides with the external momenta and the metric tensor g^''. On the left-hand 
side the resulting scalar products between the loop momentum k^ and the external 
momenta are rewritten in terms of inverse propagators, as for example 

2p-k = [-(/c-p)2 + m^] - [-fc^ + m?] + (p2 + mj-m^). (3.3) 

The first two terms of the right-hand side above cancel propagators, whereas the 
last term does not involve the loop momentum anymore. The remaining step is to 
solve for the form-factors by inverting the matrix which one obtains on the right- 
hand side of equation (|3.2I) . 

As an example we consider the two-point function: Contraction with p^ or p^j,Pu 
and 5^j. yields 

P^Bi = -- {{ml -ml- p2) Bo - Ao(mi) + ^0(^2)) , 



p' 



1 V ^21 A / -\{ml-ml-p')B^-\AQ{m2) 
P^ D ) { B22 J \ mlBo-Ao{m2) 



(3.4) 
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Figure 3 An example for an irreducible scalar product in the numerator: 
The scalar product 2p\k2 cannot be expressed in terms of inverse propagators. 



Solving for the form factors we obtain 



Bi = -—^{{ml-ml-p^)Bo-AQ{mi) + AQ{m2)) 

Zp 



B21 



1 I D , ^ 2 2n 2 D-2 ^ , ^ 
(D - l)p2 \^^"^^ -m^-p )Bi~ m^BQ —AQ[m2) 

1 



22 = 2{D~l) ^^""2 -ml- p^)Bi + 2mlBo - Ao{m2)) 



(3.5) 



Due to the matrix inversion in the last step determinants usually appear in the 
denominator of the final expression. For a three-point function we would encounter 
the Gram determinant 



A, 



PI Pi ■ P2 

Pi ■ P2 pI 



(3.6) 



One drawback of this algorithm is closely related to these determinants : In a phase 
space region where pi becomes coUinear to p2 , the Gram determinant will tend to 
zero, and the form factors will take large values, with possible large cancellations 
among them. This makes it difficult to set up a stable numerical program for au- 
tomated evaluation of tensor loop integrals. 

The Passarino-Veltman algorithm is based on the observation, that for one-loop 
integrals a scalar product of the loop momentum with an external momentum can 
be expressed as a combination of inverse propagators. This property does no longer 
hold if one goes to two or more loops. Fig. Q shows a two- loop diagram, for which 
the scalar product of a loop momentum with an external momentum cannot be 
expressed in terms of inverse propagators. 



3.2 General reduction method. Let us now consider a general multi-loop 
tensor integral and assume that we follow the steps in sect. 12.41 as we did for scalar 
integrals. After the change of variables for the diagonalisation of the quadratic 
form, we have a polynomial in the Feynman or Schwinger parameters and the loop 
momentum k in the numerator. Integrals with an odd power of the loop momentum 
in the numerator vanish by symmetry, while integrals with an even power of the 
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loop momentum can be related by Lorentz invariance to scalar integrals: 

'^''H^k'^fie) = -1-g^-^ f ^J^e)fik% (3.7) 



/ 



i^D/2 -ly J JJ^ I i^D/2 



k^^k^kPk^f{k') = ^,^ , ,^ (g^V" + g^V" + g^V") 

{-efl{k% 



i7r^/2 *' ' ' D(D + 2 

d°k 



The generalisation to arbitrary higher tensor structures is obvious. In the master 
formula eq. (|2.fi3(l we already observed that a factor (— fc^) in the numerator is 
equivalent (apart from prefactors) to a shift in the dimension D ^ D + 2. Let us 
introduce an operator D^, which shifts the dimension: 

, f d^k , ^. f d^^+'^^k , ^. , , 

°"y^^/(^') = jj^^u^wj^fik')- (3.8) 

Then we have for example 

j^.^'k-m = -y^^li^fO^')- (3-9) 

In addition, shifting the loop momentum like in k' = k — xp introduces for tensor 
integrals the (Feynman or Schwinger) parameters Xj in the numerator. For the 
tensor reduction it is convenient to work temporarily with Schwinger parameters. 
Let us recall the formula 12.29|l for Schwinger parameters: 

oo 

(Z^ - ^^ J dx x"-' eM^P)- (3.10) 



A Schwinger parameter x in the numerator is equivalent to raising the power of the 
original propagator by one unit: v ^ v + I. It is convenient to denote by i"*" the 
operator, which raises the power of propagator i by one. 



1 1 



LXJ 

: r / dXi x'^^~'^Xiexp{xiPi). (3.11) 





Therefore we can consider an integral, where a Schwinger parameter occurs in the 
numerator as a scalar integral, where the corresponding propagator is raised to a 
higher power. As a consequence, by using an intermediate Schwinger parametrisa- 
tion, we can express all tensor integrals in terms of scalar integrals [37,38]. The 
price we have to pay is that these scalar integrals involve higher powers of the 
propagators and/or have shifted dimensions. Each integral can be specified by its 
topology, its value for the dimension D and a set of indices, denoting the powers of 
the propagators. In general the number of different integrals is quite large. 

4 Scalar multi-loop integrals 

In the previous section we saw that all tensor integrals can be reduced to scalar 
integrals, where the propagators are raised to some power and the dimension is 
shifted. In sect. 12.41 we discussed these scalar integrals and gave in eq. H2.63|l a 
master formula which is applicable, once we have completed the square and shifted 
the loop momentum, such that the integrand depends only on fc^. This raises the 
question if the result of completing the square and shifting the loop momentum. 
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which has to be done iteratively for each loop, can be read off directly from the 
underlying Feynman graph. This is indeed the case and in this section we will give 
the appropriate formulae. 

4.1 Graph polynomials. We consider a scalar Z-loop integral /g in Z) = 
2to — 2e dimensions with n propagators, corresponding to the graph G: 

The momenta Qj of the propagators are linear combinations of the external and the 
loop momenta. The propagators can we raised to a power Vj . Introducing Feynman 
parameters and performing the momentum integrations as in sect. 12.41 one arrives 
at the following Feynman parameter integral [39]: 



(4.2) 



Here, v — X)7"=i ^i- The functions U and J-^ can be straightforwardly derived from 
the topology of the corresponding Feynman graph G. Cutting I lines of a given 
connected Z-loop graph such that it becomes a connected tree graph T defines a 
chord C{T,G) as being the set of lines not belonging to this tree. The Feynman 
parameters associated with each chord define a monomial of degree I. The set of all 
such trees (or 1-trees) is denoted by 7i. The 1-trees T £ Ti define U as being the 
sum over all monomials corresponding to the chords C{T, G). Cutting one more line 
of a 1-tree leads to two disconnected trees (Ti,T2), or a 2-tree. T2 is the set of all 
such pairs. The corresponding chords define monomials of degree Z + 1. Each 2-tree 
of a graph corresponds to a cut defined by cutting the lines which connected the 2 
now disconnected trees in the original graph. The square of the sum of momenta 
through the cut lines of one of the two disconnected trees Ti or T2 defines a Lorentz 
invariant 




ST ^ \ y^ P, . (4.3) 



The function J-'q is the sum over all such monomials times minus the corresponding 
invariant. The function J^ is then given by J-q plus an additional piece involving 
the internal masses rrij . In summary, the functions U and JF arc obtained from the 
graph as follows: 



U = 



E[ n ■ 

TeTi jec(T,G) 



iTi,T2)eT2 jeC{Ti.G) 

n 

T = Ta+U^Xjin'^^ . (4.4) 
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Let us consider as an example the following two-loop graph: 
Pi 

1 




For simplicity we assume that all internal propagators are massless. Then the 
functions U and J-' read: 

U = a:i5a;23 +a;i5X46 + a;23a:46, 

T = {X1X3X4 + X5X2Xe + XiX5X234:6) {~Pl) 

+ (xexsxr-, + x^xyXi + a;4a;6a;i235) (-P2) 

-I- {x2X/^x^ + xzxxxfi + 0:2X32:1456) (-P3) • (4.6) 

Here we used the notation that xij,,, — xi ^ Xj + .... In general, Z-/ is a positive 
semi-definite function. Its vanishing is related to the UV sub-divergences of the 
graph. Overall UV divergences, if present, will always be contained in the prefactor 
Y(y — ID /2). In the Euclidean region, T is also a positive semi-definite function of 
the Feynman parameters Xj . 

4.2 Feynman integrals and differential forms. Note that U, Tq and T 
are homogeneous functions in the Feynman parameters Xj with degrees 

deg U^l, AQ'gT = deg Tq = 1^\. (4.7) 

The function 

/(-I'-'-") = ln<^ I :p.-,ni. (4-8) 

is then a homogeneous function of degree 

deg/ = ~n. (4.9) 

It follows that 

/ " o \ 

+x:.,-j/ ^ 0. (4.10) 

If we now define the differential form 

n 

w = ^(-1)^""^ a:^ dxi A ... A dxj A ... A dx„, (4.11) 

where the hat indicates that the corresponding term is omitted, then /lj is closed: 

rf(/^) = 0. (4.12) 

If juj vanishes on any hyper-surface Xj = 0, then we can write the Feynman integral 
as [40] 

Ic ^ (e-.-)'£(;^i^/M (4.13) 
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where M is any hyper-surface covering the sohd angle Xj > 0. 

5 One-loop integrals 

The simplest, but most important loop integrals are the one-loop integrals. We 
have a good understanding of these integrals and I will present the main results 
in this section. For one-loop tensor integrals we can use the Passarino-Veltman 
method, which reduces any tensor integral to a combination of scalar integrals. 
Note that this method does not shift the dimension, nor does it raise the powers 
of the propagators. An important result is that any scalar integral with more 
than four external legs can be reduced to scalar integrals with no more than 4 
external legs. Therefore the set of basic one-loop integrals is rather limited. I will 
discuss this reduction in sect. 15.11 A second important result is that all one- loop 
integrals can be expressed through the logarithm and the dilogarithm. No other 
transcendental functions occur. I will discuss the appearance of the dilogarithm in 
sect. 15.21 and 15.31 Although we discussed reduction methods for tensor integrals 
already in sect.|3 the methods presented there are rather general and need not be 
the most efficient ones. Any tensor integral I^u... is usually contracted into a tensor 
jfiv... independent of the loop integration. If we have additional information on 
jA"'-- niore efficient algorithms for the tensor reduction of /^i/... can be derived. An 
example for one-loop integrals is discussed in sect. 15.41 

5.1 Reduction to integrals with no more than four external legs. In 

this section we discuss the reduction of scalar integrals with more than four external 
legs to a basic set of scalar one-, two-, three- and four-point functions. It is a long 
known fact, that higher point scalar integrals can be expressed in terms of this basic 
set [41,42], however the practical implementation within dimensional regularisation 
was only worked out recently [43-47]. The one-loop n-point functions with n > 5 
are always UV-finite, but they may have IR-divergences. Let us first assume that 
there are no IR-divergences. Then the integral is finite and can be performed in 
four dimensions. In a space of four dimensions we can have no more than four 
linearly independent vectors, therefore it comes to no surprise that in an one-loop 
integral with five or more propagators, one propagator can be expressed through 
the remaining ones. This is the basic idea for the reduction of the higher point 
scalar integrals. With slight modifications it can be generalised to dimensional 
regularisation. I will discuss the method for massless one-loop integrals 

With the notation 

i 

Qi = ^Pj (5.2) 

one can associate two matrices S and G to the integral in eq. IjS.lfl . The entries of 
the n X n kincmatical matrix S are given by 

Sij = {qi - Qjf , (5.3) 

and the entries of the (n — 1) x (n — 1) Gram matrix are defined by 

G^j = 2q,qj. (5.4) 
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For the reduction one distinguishes three different cases: Scalar pentagons (i.e. 
scalar five-point functions), scalar hexagons (scalar six-point functions) and scalar 
integrals with more than six propagators. 

Let us start with the pentagon. A five-point function in D — 4 — 2e dimensions 
can be expressed as a sum of four-point functions, where one propagator is removed, 
plus a five-point function in 6 — 2e dimensions [43]. Since the (6 — 2£)-dimensional 
pentagon is finite and comes with an extra factor of e in front, it does not contribute 
at 0{e'^). In detail we have 

5 5 

/5 = -2£B/5«-2e-^6,/f = -^Mf -KO(e), (5.5) 

1=1 j=i 

where /s"^^ denotes the (6 — 2e)-dimcnsional pentagon and l!^ denotes the four- 
point function, which is obtained from the pentagon by removing propagator i. The 
coefficients B and bi are obtained from the kinematical matrix Sij as follows: 

^' = E(0,,. B = J2b,. (5.6) 

J i 

The six-point function can be expressed as a sum of five-point functions [44] without 
any correction of 0{£) 

h = -JZb4^- (5-7) 

2=1 

The coefficients bi are again related to the kinematical matrix S^ : 

^^ - E(^"'),r (5-8) 

j 

For the seven-point function and beyond we can again express the n-point function 
as a sum over (n — l)-point functions [47]: 

n 
j=l 

In contrast to eq. H5.7|) . the decomposition in eq. (|5.9|) is no longer unique. A 
possible set of coefficients r^ can be obtained from the singular value decomposition 
of the Gram matrix 

4 

G,j = Y.^^^'^^{y^)kr (5-10) 



' kj 

fe=l 



as follows [48] 



^^ n — 1 _.n — 1 

5.2 The three-point function. As an example for the appearance of the 
dilogarithm let us discuss the one-loop three-point function with no internal masses 
and the kinematical configuration pf ^ 0, P2 ^ and p§ = (pi + P2)^ ^ 0. 

^" " '''^"^'' J z^^/2 (-p)(-(fc - pim-ik - PI - P2)2) • ^^-^^^ 
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The integral is finite and can be evaluated in four dimensions 

1 

-a;2p§ - y2p2 _^ 3.^(^2 _ p2 _ ^2) _ ^^2 _^ ^(^2 _ p2) 



Co - / dx / C?y ^2„2 _ „2„2 , ^„C„2 _ „2 _ „2^ _ ^„2 , „(^2 _ ^2\ + ^(^)- 





(5.13) 

1 follow here closely the original work of 't Hooft and Veltman [49]. We make the 
change of variables y' ~ y — ax and choose a as a root of the equation 

-a^pl + a{pl-pl-pl)-pl = 0. (5.14) 

With this choice we eliminate the quadratic term in x. We then perform the x- 
integration and we end up with three integrals of the form 
1 

— [In (ay2 + fey + c) - In (ay"^ + byo + c)] (5.15) 

y-yo 



Factorising the arguments of the logarithms, these integrals are reduced to the type 

1 
R ^ /'^^[in(y-yi)-ln(yo-2/i)]. (5.16) 

J y -yo 



This integral is expressed in terms of a new function, the dilogarithm, as follows: 

R = uJ^^)-uJj^^^^), (5.17) 

provided —yi and l/(yo ^ yi) have imaginary part of opposite sign, otherwise 
additional logarithms occur. 

5.3 The dilogarithm. The dilogarithm is defined by 

Li.(.) ^ -/rftMl_f^._|,,Ml_l) (5,3) 



If we take the main branch of the logarithm with a cut along the negative real axis, 
then the dilogarithm has a cut along the positive real axis, starting at the point 
X = I. For |x| < 1 the dilogarithm has the power series expansion 



n=l 

Some important numerical values are 

Li2(0) = 0, Li2(l) = ^, Li2(-1) = -^, Li2Q)=^-i(ln2)^ 

The dilogarithm with argument x can be related to the dilogarithms with argument 
(1 — a;) or l/x: 

Lhix) = -Li2(l-x) + -7r2-ln(.T)ln(l-a;), 

6 

Lh{x) - -Li2 Q) - ^TT^ - i (ln(-x))^ (5.20) 
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;-term relation: 








-f + Ll2 1 

1 — X J 


fxy- 
[ 1- 


-y J 


1 In^ 


' 2 \l-yj- 








(5.21) 



Li2{xy) = Li2(x) + Li2(j/) + Li2 



5.4 Spinor techniques. Tlie reduction metliods for tensor integrals discussed 
in sect. 131 are rather general and independent of the tensor structure into which the 
tensor integral is contracted. By taking into account information from this external 
tensor structure, more efficient reduction algorithms can be derived [50-55]. I will 
discuss as an example a method for one-loop integrals with massless propagators. 
The method is most conveniently explained within the FD-scheme of dimensional 
regularisation. In this scheme we can assume without loss of generality that the 
tensor structure J^i...^^^ is given by 

Jiii...p.r = (oi- |7miI^i-)---K- iTpJ^r-), (5.22) 

where (oi — | and \bj—) are Weyl spinors of definite helicity. Therefore we consider 
tensor integrals of the form 



I 



fc2(fc-pi)2...(fc-pi - ...Pn-l) 



2' 



(5.23) 



where fc['^> denotes the projection of the D dimensional vector A:^ onto the four- 
dimensional subspace. The quantity (a — |7^| 6—) is a vector in a complex vector- 
space of dimension 4 and can therefore be expressed as a linear combination of four 
basis vectors. 

The first step for the construction of the reduction algorithm based on spinor 
methods is to associate to each n-point loop integral a pair of two light-like momenta 
li and ^2, which are linear combinations of two external momenta pi and pj of the 
loop integral under consideration [53]. Obviously, this construction only makes 
sense for three-point integrals and beyond, as for two-point integrals there is only 
one independent external momentum. We write 

h^-. (pi-aipj), h = -, {~'a2Pi+Pj), (5.24) 

where ai and a2 are two constants, which can be determined from pi and pj. 

In the second step we use h and I2 to write (a — |7^| b—) as a linear combination 
of the basis vectors 

{h-\j^\k-), (?2_|7^|/2-), (z, _ I^^J /2-) , {l2-\l,\li-). (5.25) 

The contraction of k'f^s with the first or second basis vector leads to 

'^kih = ['2piki ~ ai2pjki] , 2fc//2 = :; [-a2'2piki + 2pjki] , 

i — aia2 t — aia2 

(5.26) 

and therefore reduces immediately the rank of the tensor integral. Repeating this 
procedure we end up with integrals, where the numerator is given by products of 

(Zi- 1^(4)1 /2-) and (?2- 1^(4)1^1-), (5.27) 

plus additional reduced integrals. Therefore the tensor integral is now in a standard 
form. In the next step one shows that any product of factors as in eq. (|5.27|l can 
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V 
1 2^ 

3 



Figure 4 The triangle trick: Integration by parts eliminates the propagators 
1, 2, 3 or 4. 



be reduced. This reduces all tensor integrals to rank 1 integrals. Finally, the rank 
1 integrals are reduced to scalar integrals. 

6 Advanced methods 

In sects. |5] to 0] we discussed general /-loop integrals and showed that all ten- 
sor integrals can be reduced to scalar integrals and that the integration over the 
I independent loop momenta can be performed at the expense of introducing ad- 
ditional parameter integrals. The integrand for the Feynman parameter integral 
can be read off directly from the underlying Feynman graph. In this section I will 
present methods to compute these Feynman parameter integrals. As the tensor 
reduction usually introduces a huge number scalar integrals, which differ by the 
powers of the propagators or the by the dimension, I start in sect. 16. II and 16. 21 with 
techniques, which reduce the unknown integrals to a smaller set, called "master 
integrals". In sects. IH7^ to 16. 51 1 will present techniques for the analytic calculation 
of Feynman parameter integrals. Finally I review in sect. 16.61 a method for the 
numerical computation of the coefhcients of the Laurent expansion. 

6.1 Integration by parts. Integration-by-part identities [56] are based on 
the fact, that the integral of a total derivative is zero: 

-v^f{k,p,) = 0. (6.1) 



Here, k is the loop momentum, the p^'s are external momenta with respect to this 
loop integration and v can either be a loop momentum or an external momentum. 
Working out the derivative yields a relation among several scalar integrals. 

As an example we look at the triangle trick [57]: Assume that we have inside 
a multi-loop integral the building block shown in fig. ^ with massless propagators 
and p^ = 0: 

T{vi,l'2,V'i,l'A,Vz) = 



{-kir{-{k,-pYY^ 

d^ko 1 



*7r^/2 {-kir {-{k2~p?r {-{k2 - kiYr ■ 

Then by choosing k = k2 and v = k2 in eq. H6.1|l we obtain 



(6.2) 



The choice k — k2 ~ p and v — k2 ~ p yields a second relation: 

[{D -V2- 21^3 - Vo) - i^22+3- - 1/55+ (3- - 4")] T{vx, V2, v^., 1/4, t^s) = 



(6.3) 



(6.4) 
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If we assume that vi and 1/2 are positive integers we can eliminate by repeated 
use of eq. H6.3|l either propagator 1 or 2. Note that if propagator 2 is eliminated, 
propagators 3 and 5 form a rather trivial one-loop insertion. Similar considerations 
apply to propagators 3 and 4 by using eq. H6.4(l . 

6.2 Reduction to master integrals. In the example discussed above we 
could use the triangle trick to simplify the integrand. In general, integration-by- 
part identities lead to relations among various scalar integrals. Additional relations 
among scalar integrals are obtained from the invariance of scalar integrals under 
Lorentz transformations [58]. A scalar integral is evidently invariant under an 
infinitesimal Lorentz transformation, parametrised as 

p" ^p'' + Spf" ^p'' + Sei^p" with <5e;; = -St^ . (6.5) 

This implies that 

Pl^ Pl^ ^ ■■■+Pn^ Pn^ ) I{Pl,---,Pn) = , (6.6) 

where I{pi, ■■■,Pn) is a scalar integral with the external momenta Pi. Again, work- 
ing out the derivatives yields a relation among several scalar integrals. 
Each relation is linear in the scalar integrals and in principle one could use Gauss 
elimination to reduce the set of scalar integrals to a small set of master integrals. 
In practice this approach is inefficient. An efficient algorithm has been given by La- 
porta [59]. The starting point is to introduce an order relation for scalar integrals. 
This can be done in several ways, a possible choice is to order the topologies ffist: 
A scalar integral corresponding to a topology Ti is considered to be "smaller" than 
an integral with topology T2, if Ti can be obtained from T2 by pinching of some 
propagators. Within each topology, the scalar integrals can be ordered according 
to the powers of the propagators and the dimension of space-time. Laporta's algo- 
rithm is based on the fact that starting from a specific topology, integration-by-part 
and Lorentz- invariance relations only generate relations involving this topology and 
"smaller" ones. To avoid to substitute a specific identity into a large number of 
other identities, one starts from the "smallest" topology and generates all relevant 
relations for this topology. Inside this class, integrals with higher powers of the 
propagators or higher dimensions are then expressed in terms of a few master in- 
tegrals. These manipulations involve only a small subset of the complete system 
of integrals and relations and can therefore be done efficiently. Once this topology 
is completed, one moves on to the next topology, until all topologies have been 
considered. All integrals, which cannot be eliminated by this procedure are called 
master integrals. These master integrals have to be evaluated explicitly. 

6.3 Mellin-Barnes representation. In sect. 0] we saw that the Feynman 
parameter integrals depend on two graph polynomials U and J-', which are homo- 
geneous functions of the Feynman parameters. In this section we will continue the 
discussion how these integrals can be performed and exchanged against a (multiple) 
sum over residues. The case, where the two polynomials are absent is particular 
simple: 

n 

nr(^.) 
r(j/i + ... + !/„)' 
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With the help of the Melhn-Barnes transformation we now reduce the general case 
to eq. H6.7|l . The Mellin-Barnes transformation reads 

ioo ioo 

{Ai +A2 + ... + An^ ^ ] ]^_^ [ dai... [ da^-i (6.8) 



— 'ioc 



(Ti A'^ri-1 A-CTi-...-(Tn-l-C 



xr(-ai)...r(-a„_i)r(ai + ... + an-i + c) A^' ...A^^SiA- 

Each contour is such that the poles of r(— tr) are to the right and the poles of r(cr+c) 
are to the left. This transformation can be used to convert the sum of monomials of 
the polynomials ZY and J^ into a product, such that all Feynman parameter integrals 
are of the form of eq. (|6.7|) . Therefore we exchange the Feynman parameter integrals 
against multiple complex contour integrals. As this transformation converts sums 
into products it is the "inverse" of Feynman parametrisation. The contour integrals 
are then performed by closing the contour at infinity and summing up all residues 
which lie inside the contour. Here it is useful to know the residues of the Gamma 
function: 

(-1)" (-1)" 
res (r(cr + a),cr = —a— n) — j — , res (r(— ct + a),a = a + n) = j — . 



(6.9) 

Therefore we obtain (multiple) sum over residues. Techniques to manipulate these 
sums are discussed in the next section. In particular simple cases the contour 
integrals can be performed in closed form with the help of two lemmas of Barnes. 
Barnes first lemma states that 



-^ / d<7T(a + a)T(b + a)T(c-a)T(d~(T) 
Im J 



V{a + c)T{a + d)T{b + c)V{b + d) 
2711 J "'""^"' ' -y-v" ' -y-v- -y-v- -/ T{a + b + c + d) ' 

—ioo 

(6.10) 

if none of the poles of r(a + cr)r(fo+cr) coincides with the ones from V{c—a)T{d—a). 
Barnes second lemma reads 

ioo 

1 f J r(a + cr)r(6 + cr)r(c + cr)r(d - cr)r(e - cr) 
da- 



(6.11) 



2m J T{a + b + c + d + e + a) 

— ioc 

Tja + d)r{b + d)r{c + d)r{a + e)r(6 + e)r(c + e) 
~ r{a + b + d + e)r{a + c + d + e)r{b + c + d + e) ' 



6.4 Expansion in a small parameter for specific classes of transcen- 
dental functions. From the Mellin-Barnes representation we obtain multiple sums 
involving Gamma functions. The small parameter e appears in the arguments of 
the Gamma functions. We are interested in the Laurent expansion in e. In this 
section we discuss techniques to manipulate multiple sums. As an example we con- 
sider again the scalar integral already encountered in eq. (|2.62|l . but this time in 
dimension D = 2m — 2e and each propagator Pj is raised to the power Vj . A short 
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calculation leads to 

i2m-2et. Ill 1 

j^m-e (_fc2)., (_^2)., (_fc2).3 I ^123^ T{vi)T{v2) ^ ' 

r(m - £ - ;/i)r(m - £ - t/23) v^ r(n + t^2)r(n - m + £ + t/123) , _ .„ 
"^ r(2m - 2£ - i.123) ^0 r(n+l)r(n + ^.23) ^ ""^ ' 

where x = S12/S123, 1^23 = ^"2 + 1^3 and 2^123 = i^i + z/2 + J^3- The infinite sum in 
the last line of (|6.12l) is a hyper-geometric function, where the small parameter e 
occurs in the Gamma functions. More complicated loop integrals yield additional 
classes of infinite sums. One often encounters the following types of infinite sums: 
Type A: 



E 

1=0 



r(i + ai) r(z + afc) 



x' (6.13) 



Up to prefactors the hyper-geometric functions j+iFj fall into this class. The 
example discussed above is also contained in this class. 
Type B: 



00 OD 



T{i + ai) r(i-hafe)r(i + 6i) T{j + bi) 



~^^4±1±^},J4±1±^ xW (6.14) 

An example for a function of this type is given by the first Appell function Fi . 
Type C: 

Here, an example is given by the Kampe de Feriet function Si. 
Type D: 

hhS ^ )n^ + <y■n^+a',)T{J+b',)■■■T{j+h[) 



^r(. + j + ci)-r(. + , + c)"^ (^-^^^ 

An example for a function of this type is the second Appell function F2 . 
Note that in these examples there are always as many Gamma functions in the 
numerator as in the denominator. We assume that all a„, a^, &„, 6^, Cn and c^ are 
of the form "integer -I- const • £" . The generalisation towards the form "rational 
number + const • £" is briefly discussed in sect. 16.4.31 The task is now to expand 
these functions systematically into a Laurent series in e. 

6.4.1 Nested sums. In this section I review the underlying mathematical struc- 
ture for the systematic expansion of the functions in H6.13|l - (|6.16|l . I discuss proper- 
ties of particular forms of nested sums, which are called Z-sums and show that they 
form a Hopf algebra [60-62] . This Hopf algebra admits as additional structures a 
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conjugation and a convolution product. Z-sums are defined by 

Z{n;mi,...,mk;xi,...,Xk) = V -^•••-^- (6-17) 

— *1 ik 
n>ii>i2>...>ik>0 

k is called the depth of the Z-sum and w = mi + ... + mk is called the weight. If 
the sums go to Infinity (n = oo) the Z-sums are multiple polylogarithms [63]: 

Zioo;mi,...,mk;xi,...,Xk) = Li„j,...,„^(a:;i, ..., Xfe). (6.18) 

For xi — ... — Xk — 1 the definition reduces to the Euler-Zagier sums [64,65]: 

Z(n;mi,...,mfe;l,...,l) = Zmi,...,mfc(«)- (6.19) 

For n = oo and xi ^ ... = Xk = I the sum is a multiple (^- value [66]: 

Z(oo;mi,...,mfc;l,...,l) == Cmi,...,mk- (6.20) 

The usefulness of the Z-sums lies in the fact, that they interpolate between multiple 
polylogarithms and Euler-Zagier sums. 

In addition to Z-sums, it is sometimes useful to introduce as well S'-sums. 
5-sums are defined by 

S{n;mi,...,mk;xi,...,Xk) = V -^•••-^- (6-21) 

— 1l Ik 
n>ii>i2>...>ifc>l 

The S'-sums reduce for xi = ... = .Xfe == 1 (and positive mi) to harmonic sums 
[67-69]: 

S'(n;mi,...,mfe;l,...,l) = Srm.,...,mk{'>^)- (6.22) 

The S-sums are closely related to the Z-sums, the difference being the upper sum- 
mation boundary for the nested sums: (i — 1) for Z-sums, i for S'-sums. The 
introduction of S'-sums is redundant, since S-sums can be expressed in terms of 
Z-sums and vice versa. It is however convenient to introduce both Z-sunis and 
S-sums, since some properties are more naturally expressed in terms of Z-sums 
while others are more naturally expressed in terms of S-sums. An algorithm for 
the conversion from Z-sums to S-sums and vice versa can be found in [60]. 

The Z-sums form an algebra. The unit clement in the algebra is given by the 
empty sum 

e = Z{n). (6.23) 

The empty sum Z{n) equals 1 for non-negative integer n. Before I discuss the mul- 
tiplication rule, let me note that the basic building blocks of Z-sums are expressions 
of the form 

re" 
-^, (6.24) 

which will be called "letters". For fixed n, one can multiply two letters with the 
same n: 

x^x^^ {x,x,r (g25) 

e.g. the Xj's are multiplied and the degrees are added. Let us call the set of all 
letters the alphabet A. As a short-hand notation I will in the following denote a 
letter just by Xj — x'^/n™^ . A word is an ordered sequence of letters, e.g. 

W = Xi,X2,...,Xk. (6.26) 



30 



«2 i 

















Stefan ^Veinzierl 




«2 1 




»2 1 




«2 1 


• • • • 










• • • 






• 


• • • • 




• 






• • 






• 


• • • • 


= 


• • 


+ 




• 


+ 




• 


• • • • 




• • • 












• 



n 



H 



n 



«i 



Figure 5 Sketch of the proof for the multiplication of Z-sums. The sum over 
the square is replaced by the sum over the three regions on the r.h.s. 



The word of length zero is denoted by e. The Z-sums defined in (|6.17|) are therefore 
completely specified by the upper summation limit n and a word W. A quasi-shuffle 
algebra A on the vector-space of words is defined by [70] 



eoW 



Woe = W, 

Xi, {Wi o (^2, W2)) + X2, {{X,,Wi) o W2) 

+ {Xi ■ X2), {Wi o W2). (6.27) 

Note that "•" denotes multiplication of letters as defined in eq. (|6.25() . whereas "o" 
denotes the product in the algebra A^ recursively defined in eq. lfi.27|l . This defines 
a quasi-shuffle product for Z-sums. It is also called a mixable shuffle product [71] 
or stuffle product [66] . The recursive definition in H6.27|l translates for Z-sums into 






1 



/, -nil ^ni2 
1*1 



il 



■E 

12 = 1 *2 



~^^'n 



i—\ 



T^n 



Xh - l)Z,n' ,n'{i2 - 1) 



,,(i- l)Z„,.^,...,^;(i- 1). 



(6.28) 



The proof that Z-sums obey the quasi-shuffle algebra is sketched in Fig. The 
outermost sums of the Z-sums on the l.h.s of H().28|l are split into the three regions 
indicated in Fig. [S] A simple example for the multiplication of two Z-sums is 

Z{n;mi;xi)Z{n;mi;x2) — (6.29) 

Z(n; TTii, TO2; Xi, X2) + Z(n; m2, mi, X2, Xi) + Z{n; mi + 1122; a;ia;2). 

The quasi-shuffle algebra A is isomorphic to the free polynomial algebra on 
the Lyndon words. If one introduces a lexicographic ordering on the letters of the 
alphabet A, a Lyndon word is defined by the property 

W <V (6.30) 

for any sub- words U and V such that W = U,V. Here U, V means just concatena- 
tion of U and V. 

The Z-sums form actual a Hopf algebra. It is convenient to phrase the coalgebra 
structure in terms of rooted trees. Z-sums can be represented as rooted trees 
without any side-branchings. As a concrete example the pictorial representation of 
a sum of depth three reads: 
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n ii —1 22 — 1 



Z{n;mi,m2,m3;xi,X2,X3) = X! X] 51 ~^ 



1 12 = 1 13 = 1 



Xi 
X2 
X3 



(6.31) 




The outermost sum corresponds to the root. By convention, the root is always 
drawn on the top. Trees with side-branchings are given by nested sums with more 
than one sub-sum, for example: 



^T^^i'i'~'^'^'^2,X2)Z{i-l;m3;x3) = / \ (6.32) 

i=l 

Of course, due to the muhiphcation formula, trees with side-branchings can always 
be reduced to trees without any side-branchings. The coalgebra structure is now 
formulated in terms of rooted trees. 1 first introduce some notation how to ma- 
nipulate rooted trees, following the notation of Kreimer and Connes [72,73]. An 
elementary cut of a rooted tree is a cut at a single chosen edge. An admissible cut 
is any assignment of elementary cuts to a rooted tree such that any path from any 
vertex of the tree to the root has at most one elementary cut. An admissible cut 
maps a tree i to a monomial in trees ii o ... o tk+i- Note that precisely one of these 
subtrees tj will contain the root of t. Denote this distinguished tree by R'-^ (t), and 
the monomial delivered by the k other factors by P'~^{t). The counit e is given by 

e(e) = 1, 

e{t) = 0, t^e. (6.33) 

The coproduct A is defined by the equations 

A(e) = e(g)e, 

A{t) = e(E)t + t(E)e+ ^ P^{t)®R^{t), 

adni. cuts O of t 

A(iio...oifc) = A(ii)(o®o)...(o(g)o)A(ifc). (6.34) 

The antipode S is given by 

5(e) = e, 

S{t) = t J2 S{P^{t))oR^{t), 

adm. cuts O of t 

S{tio...otk) = S{ti)o ...oS{tk). (6.35) 

Since the multiplication in the algebra is commutative the antipode satisfies 

S^ = id. (6.36) 

Let me give some examples for the coproduct and the antipode for Z-sums: 
AZ{n;mi;xi) = e (^ Z{n;mi;xi) + Z{n;mi;xi) (^ e, 
AZ{n;mi,m2;xi,X2) — e ^ Z{n;mi,m2;xi,X2) + Z{n;mi,m2;xi,X2) ^ e 

+Z{n;m2;x2) ^ Z{n;mi;xi), (6.37) 

SZ{n;mi;xi) = —Z{n;mi;xi), 
SZ{n;mi,m2;xi,X2) = Z{n;m2,mi;x2,xi) + Z{n;mi + m2;xiX2)- 
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The Hopf algebra of nested sums has additional structures if we allow expressions 
of the form 

— 2_2'(n;TOi, ...,TOfc;a;i, ...,Xk), (6.38) 

e.g. Z-sums multiplied by a letter. Then the following convolution product 



n-l 



,,n—i 



E |^^(^ - 1' ■•■) 7;r:^^(" - * - 1; ■•■) (6.39) 

can again be expressed in terms of expressions of the form (|6.38|) . An example is 

5^ -Zl(^-l)-LZl(r^-^-l) = 
i=i ^ ' 



X 

n 



Z ( n - 1; 1, 1, 1; ^, -, ^ ) + Z ( n - 1; 1, 1, 1; ^, 1, - 



+Z(n-l;l, 1,1; 1,^,1 
V X 

In addition there is a conjugation, e.g. sums of the form 



X y xj \ X y ^ 

+ {x^y). (6.40) 



-E U-^y T^S{^■,...) (6.41) 



« / I 



i=l 

can also be reduced to terms of the form (|6.38|) . Although one can easily convert 
between the notations for S'-sunis and Z-sums, expressions involving a conjugation 
tend to be shorter when expressed in terms of 5-sums. The name conjugation stems 
from the following fact: To any function f{n) of an integer variable n one can define 
a conjugated function C o f{n) as the following sum 

Co fin) = J2( ^)(-l)7(*). (6.42) 

Then conjugation satisfies the following two properties: 

Col = 1, 
CoCofin) = fin). (6.43) 

An example for a sum involving a conjugation is 



-E : (-1) T^^« 



S (n- 1, 1; 1 - X, -—) ~ S in; 1, 1; 1 - x, 1) . (6.44) 

Finally there is the combination of conjugation and convolution, e.g. sums of the 
form 

ft J./ \ T in — ^ 



4=1 



E( "j(-i)V^(^--)(;rr7F^(^^~^'-) ('-^'^ 
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can also be reduced to terms of the form (|6.38|l . An example is given by 
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n-1 

E 



(-1)' S{i;l;x) S{n-i;l;y) 



(6.46) 



Sin;l;y) + il~xy 



S n:l 



(-1)" 
+^^)s(n;l;x) + il-yy 



1-i 



S n: 1; 



y\i 



1-- 



S \ n:l: 



6.4.2 Expansion of hyper- geometric functions. In this section I discuss how the 
algebraic tools introduced in the previous section can be used for the Laurent ex- 
pansion of the transcendental functions (|6.13() - (|6.16() . First I give some motivation 
for the introduction of Z-sums: The essential point is that Z-sums interpolate be- 
tween multiple polylogarithms and Euler-Zagier-sums, such that the interpolation 
is compatible with the algebra structure. On the one hand, we expect multiple 
polylogarithm to appear in the Laurent expansion of the transcendental functions 
(I6.13() - (|6.16() . a fact which is confirmed a posteriori. Therefore it is important that 
multiple polylogarithms are contained in the class of Z-sums. On the other the 
expansion parameter e occurs in the functions H6.13|l - H6.16|l inside the arguments of 
Gamma- functions. The basic formula for the expansion of Gamma functions reads 



r(n + e) = r(l + e)T{n) [l + eZi{n - 1) + e^Zii{n - 1) 
+e^Zni{n - 1) + ... + e"-iZn...i(n - 1)] , 



(6.47) 



containing Euler-Zagier sums for finite n. As a simple example I discuss the expan- 
sion of 



E 



r(z + ai+ tie)T{i + Qz + tag) 
T{i + l)T{i + as + tse) 



(6.48) 



into a Laurent series in e. Here oi, 02 and 03 are assumed to be integers. Up 
to prefactors the expression in H6.48|l is a hyper-geometric function 2Fi- Using 
r(a; -|- 1) = xT{x), partial fractioning and an adjustment of the sununation index 
one can transform (16.481) into terms of the form 



E 



T{i + tie)T{i + t2e) x' 



r(?;)r(i + tse) i™' 

where m is an integer. Now using 16.4711 one obtains 

(1 + etiZi{i - 1) + ...) (1 + et2Zi{i - 1) + ...) x' 



r(i + £)E' 



{l + et3Zi{i-l) + ...) 



(6.49) 



(6.50) 



Inverting the power series in the denominator and truncating in e one obtains in 
each order in e terms of the form 



E' 

4=1 



,Ji - l)Z,„/...„i/(i - l)2'„i/'...„j//(i - 1) 



(6.51) 
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Using the quasi-shuffle product for Z-sunis the three Euler-Zagier sums can be 
reduced to single Euler-Zagier sums and one finally arrives at terms of the form 

oo ^ 

J2^Z^l-rnA^-^), (6.52) 

i=l 

which are special cases of multiple polylogarithms, called harmonic polylogarithms 
Hm,mi....,mk{x)- This Completes the algorithm for the expansion in e for sums of the 
form Hfi.lHfl . Since the one-loop integral discussed in 1)6.12(1 is a special case of (|6.1HII , 
this algorithm also applies to the integral H6.12I) . In addition, this algorithm shows 
that in the expansion of hyper-geometric functions j+iFj{ai, ..., aj+i; 6i, ..., bj; x) 
around integer values of the parameters a^ and bi only harmonic polylogarithms 
appear in the result. 

The algorithm for the expansion of sums of type ((6.13(1 used the multiplication 
formula for Z-sums to pass from ((6.51(1 to ((6.52(1 . To expand double sums of type 
((6.14() one needs in addition the convolution product 1(6.39(1 . To expand sums of 
type 1(6.15(1 the conjugation 1(6.41(1 is needed. Finally, for sums of type 1(6.16(1 the 
combination of conjugation and convolution as in 16.45(1 is required. More details 
can be found in [60]. 

Let me come back to the example of the one-loop Feynman integral discussed 
in eq. ((2.62(1 . For i^i = 1/2 = J^s = I and m = 2 in 1(6.12(1 one obtains: 



j4-2e 



ki 1 



•2- i-kf) i-kl) i^ki) 

r(-£)r(i-e)r(i + e)(-p?23) 




(6.53) 
1(1 -a;). 



r(l-2e) 1-x 

In this special case all harmonic polylogarithms can be expressed in terms of powers 
of the standard logarithm: 

^l,...,l(l-^) = ^^(Ina;)". (6.54) 



This particular example is very simple and one recovers the well-known all-order 
result 

r(l^e)^r(l+e)(-pL3)'''^l~^- ,,,,. 

r(i - 2e) 7^ r^^' ^^■^^> 

which (for this simple example) can also be obtained by direct integration. 

6.4.3 Expansion around half-integer values. Up to now we only discussed the 
case where the Gamma functions are expanded around an integer value. The exten- 
sion to rational numbers is straightforward, if the Gamma functions always occur 
in ratios of the form 

r(n + a - £ + be) 

r(n + c-H+&) ' (^-'^^ 
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where the same rational number p/q occurs in the numerator and in the denomi- 
nator [74] . The generahsation of eq. H6.47(l reads 



r( n+l--+£ 



r(i-H + ,)r(n + i-H) 
r(i-f) 

^ /=0 k=l 




and introduces the q-th roots of unity 

rl = e^^n^^Y (6.58) 

With the help of the q-th roots of unity we may express any Z-sum Z{n] ...) as a 
combination of Z-sums Z{q ■ n; ...), where the summation goes now up to g • n. If 
the occurrence of rational numbers is not balanced as in eq. H6.56|l . the sums can 
be performed only in special cases. These include the important cases of binomial 
sums 

E( ^ ] ^S{n;mi,...,mk;xi,...,Xk), (6.59) 



and inverse binomial sums 



y^±^^. (6.60) 

n— 1 

V ^ 

6.5 Differential equations. An alternative approach to the computation of 
Feynman parameter integrals is based on differential equations [75-77] . To evaluate 
these integrals the procedure used in [58, 78, 79] consists in finding first for each 
master integral a differential equation, which this master integral has to satisfy. 
The derivative is taken with respect to an external scale, or a ratio of two scales. 
It turns out that the resulting differential equations are linear, inhomogeneous first 
order equations of the form 

^T{y) + f{y)T{y)^g{y), (6.61) 

where y is usually a ratio of two kinematical invariants and T(jj) is the master inte- 
gral under consideration. The inhomogeneous term g{y) is usually a combination of 
simpler master integrals. In general, a first order linear inhomogeneous differential 
equation is solved by first considering the corresponding homogeneous equation. 
One introduces an integrating factor 

M{y) = e'l'f'^y^'^y, (6.62) 

such that To{y) — 1/M{y) solves the homogeneous differential equation {g{y) = 0). 
This yields the general solution of the inhomogeneous equation as 



Tiy) = j^ ^J 9{y)M{y)dy + Cj , (6.63) 

where the integration constant C can be adjusted to match the boundary condi- 
tions. This method yields a master integral in the form of transcendental functions 
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(e.g. for example hyper-geometric functions), which still have to be expanded in 
the small parameter e of dimensional regularisation. Although this can be done sys- 
tematically and was discussed above, the authors of [58,78,79] followed a different 
road: Factoring out a trivial dimension-full normalisation factor, one writes down 
an ansatz for the solution of the differential equation as a Laurent expression in 
e. Each term in this Laurent series is a sum of terms, consisting of basis functions 
times some unknown (and to be determined) coefficients. This ansatz is inserted 
into the differential equation and the unknown coefficients are determined order 
by order from the differential equation. This approach will succeed if we know in 
advance the right set of basis functions. The basis functions are taken as a subset 
of multiple poly logarithms. 

6.6 Sector decomposition. In sect. 0] we presented the general formula for 
a scalar Z-loop integral as a Feynman parameter integral. We are interested in the 
Laurent expansion in e of this integral. In this section I will discuss an algorithm, 
which allows to compute the coefficients of the Laurent expansion numerically for 
a given kinematical configuration of external momenta. The major challenge such 
an algorithm has to face is the disentanglement of overlapping singularities. An 
example for an overlapping singularity is given by 

The term l/{xi -1-2:2) is an overlapping singularity. Sector decomposition [80-82] 
is a convenient tool to disentangle overlapping singularities. For the example in 
eq. (|6.64|l one splits the integration region into two sectors xi > X2 and xi < X2- 
In the first sector one rescales X2 as x'2 — 2:2 /xi, while in the second sector one 
rescales x'l = X\lx2- Binoth and Heinrich [82] gave a systematic algorithm for 
integrals of the form of eq. (|4.2|l . First the integration domain is mapped from the 
(n — l)-dimensional simplex to the (n — l)-dimensional hypercube. This is done as 
follows: First the integration domain is split into n parts, using the identity 

/ d"x = 51/ d'^xJleixi-Xj) , (6.65) 

■^0 1=1 Jo 

such that the integral Iq becomes a sum over n integrals Iq , where in each "primary 
sector" / the variable xi is the largest one. The variables are transformed in each 
primary sector / as follows: 

J <l, 

J = I, (6.66) 

By construction, xi factorises from U and J-', xi is eliminated in each Gi using 

n-l 

dxi/xi6{l-xi{l + Y,tk))^l- (6.67) 

fe=i 

This ensures that the singular behaviour leading to poles in e still comes from the 
region of small i's. By applying the sector decomposition iteratively, one finally 
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arrives at a form where all singularities are factorised explicitly in terms of factors 
of Feynman parameters like t~ ~'"^. Subtractions of the form 

dt, tj'-'^' fit,) = -- /(O) + / dt, tj'-^' [fit,) - /(O)] (6.68) 

'*£ Jo 

for each j, where limj.^o fi^j) is finite by construction, allow to extract all poles 
and lead to integrals which are finite and can be integrated numerically. 

7 Multiple polylogarithms 

In this section I discuss multiple polylogarithms. The multiple polylogarithms 
are defined by 

Li„H....,m,(a:i,...,Xfe) = ^ -^■■■-4^- C^-l) 

ll>l2>->»l->0 

They are special cases of Z-sums: 

Limi,....mfc(a;i,...,Xfe) = Z(cx);mi, ..., mfe;xi, ..., a;^). (7.2) 

The multiple polylogarithms contain as the notation already suggests as subsets the 
classical polylogarithms Li„(a;) [83], as well as Nielsen's generalised polylogarithms 

[84] 

Sn,pix) = Li„+i_i,...,i(a;, !,...,!), (7.3) 

and the harmonic polylogarithms [85] 

^nii,...,'mk\-^ ) — ^^7n-i.....7nk\'^T ^j '••7 ^ ) ' V' J 

fe-1 

Multiple polylogarithms have been studied extensively in the literature by physicists 
[19,67,78,85-88] and mathematicians [63,66,89-98]. Here I summarise the most 
important properties. Being special cases of Z-sums they obey the quasi-shuffle 
Hopf algebra for Z-sums. In addition there is a second algebra structure, derived 
from the iterated integral representation. 

7.1 Definition through iterated integrals. Multiple polylogarithms have 
been defined in this article via the sum representation l|7.1|l . In addition, they 
admit an integral representation. From this integral representation a second algebra 
structure arises, which turns out to be a shuffle Hopf algebra. To discuss this second 
Hopf algebra it is convenient to introduce for Zk ^ the following functions 

G(.„...,z.;,) . Jdi^Jdt^ydt^^ (,.5) 

J ti- Zi J t2- Z2 J tk- Zk 



In this definition one variable is redundant due to the following scaling relation: 

G{zi,--;Zk;y) = Gixzi, ...,xzk;xy) (7.6) 

If one further defines 

9{z;y) = -^, (7.7) 

y- z 
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then one has 

d 



G{zi,...,zk;y) = g{zi;y)G{z2, ...,Zk;y) (7.. 

ay 



and 



y 
G{zi,Z2,...,Zk;y) = dt g{zi;t)G{z2,...,Zk;t). (7.9) 



One can shghtly enlarge the set and define G(0, ..., 0; y) with k zeros for zi to zt to 
be 

G(0,...,0;2/) = ;^(lny)^ (7.10) 

This permits us to allow trailing zeros in the sequence (zi, ...,Zfc) by defining the 
function G with trailing zeros via H7.9(l and (|7.1()|l . To relate the multiple polylog- 
arithms to the functions G it is convenient to introduce the following short-hand 
notation: 

Gmi,...,7nk{zi,-;Zk;y) = 0(0, . .. , 0, Zi , . .. , 2fc_i , 0.. ., 0, Zfc; y) (7.11) 

irii-1 nifc— 1 

Here, all Zj for j = 1, ..., /c are assumed to be non-zero. One then finds 

Llmi,...,mfc(2;i, ..., Xfc) = (^1) Gmi,....mk \ 7 j---? il)-(7.12) 

\Xl 0:1X2 Xl...Xk J 

The inverse formula reads 

Grm,...,mk{zi,---,zk;y) = (-I)'' Li„, „J —,—,..., -^^ ) . (7.13) 

V^^l Z2 Zk ) 

Eq. (|7.12|) together with (|7.11|) and H7.5|l defines an integral representation for the 
multiple poly logarithms. To make this more explicit I first introduce some notation 
for iterated integrals 

A t„ t2 

(7.14) 



(7.15) 



A A 

( dt dt r dt, 


■I f dtn-1 ^^ ^^ f dti 


J t~ an t - ai J tn" 



an J tji-l — fln-l J ti — ai 



and the short hand notation: 




f [dt\"' dt 


A 

r dt dt dt 


J \tj t-a 




J J "^ ■'■ t\-a' 
m times 


The integral representation for Limj... 


^mj.(xi, ..., Xfe) reads then 


Limi,...,mfc(2:ij ■■■,Xk) = 


^->'/(T»r^ 




fdt \"''-^ dt 
o — o 

\t J t- 


fdt y^"-^ dt 

-— o ... o — o -—^ 

62 \t ) t^bk' 


where the 6j's are related to the Xj's 




bj = 


1 


XiX2...Xj 



(7.16) 



(7.17) 
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Figure 6 Sketch of the proof for the multiphcation of two G-functions. The 
integral over the square is replaced by two integrals over the upper and lower 
triangle. 



From the iterated integral representation H7.5(l a second algebra structure for the 
functions G{zi, ..., Zk;y) (and through H7.12|l also for the multiple polylogarithms) 
is obtained as follows: We take the Zj's as letters and call a sequence of ordered 
letters w — zi, ...,Zk a word. Then the function G{zi, ...,Zk]y) is uniquely specified 
by the word w = zi, ...,Zfc and the variable y. The neutral element e is given by 
the empty word, equivalent to 



G{-y) = 1. 
A shuffle algebra on the vector space of words is defined by 



(7.18) 



eo w 



woe 



w, 



{zi,Wl)o {Z2,W2) = Zi,{wiO {Z2,W2)) + Z2,{{zi,Wl)oW2)- (7.19) 

Note that this definition is very similar to the definition of the quasi-shuffle algebra 
(I6.27|l . except that the third term in l|6.27|l is missing. In fact, a shuffle algebra is a 
special case of a quasi-shuffle algebra, where the product of two letters is degenerate: 
Xi ■ X2 — for all letters Xi and X2 in the notation of Sect. 16.4.11 The definition 
of the shuffle product (|7.19|) translates into the following recursive definition of the 
product of two G-functions: 



G{zi, ...,zk;y) x G{zk+i, 
y 



dt 



t ~ Zl 



G(Z2,. 



Zn]y) = 

zk;t)G{zk+i, . 



1 Zri: ^J 



-I- 



dt 



t — Zk+l 



G{zi,...,Zk;t)G{zk 



+2, 



■,t) 



(7.20) 



The proof is sketched in fig. IHl 

For the discussion of the coalgebra part for the functions G(2:i , . . . , z^ ; y) we may 
proceed as in Sect. IS. 4. H and associate to any function G(zi, ..., z^; y) a rooted tree 
without side-branchings as in the following example: 



Zlf 

G{zi,Z2,Z3;y) = Z2 

Z3i 



(7.21) 



The outermost integration (involving zi) corresponds to the root. The formulae 
for the coproduct H6.34|l and the antipode H6.^^5|l apply then also to the functions 
G{zi,...,zk;y). 
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A shuffle algebra is simpler than a quasi-shuffle algebra and one finds for a 
shuffle algebra besides the recursive definitions of the product, the coproduct and 
the antipode also closed formulae for these operations. For the product one has 

G{zi,...,zk;y) G{zk+i,...,zk+i;y) = ^ G (z^(i), ...,z<^(fe+,);y) , (7.22) 

shuf fie 

where the sum is over all permutations which preserve the relative order of the 
strings zi, ...,Zk and Zfc+i, ■■■,Zk+i- This explains the name "shuffle product". For 
the coproduct one has 

k 

AG{zi,...,Zk;y) = ^G{zi,...,zf,y)®G{zj+i,...,Zk]y) (7.23) 

and for the antipode one finds 

5G(zi,...,^fc;2/) = (-l)'=G(zfc,...,zi;y). (7.24) 

The shuffle multiplication is commutative and the antipode satisfies therefore 

S^ = id. (7.25) 

From (|7.24(l this is evident. 

7.2 The antipode and integration-by-parts. Integration-by-parts has al- 
ways been a powerful tool for calculations in particle physics. By using integration- 
by-parts one may obtain an identity between various G- functions. The starting 
point is as follows: 

y 

G{zi, ..., zk-, y)^ Jdt ('^G(zi; i) j G(z2, ..., Zk;y) (7.26) 



y 
= G{zi;y)G{z2, ..., zk;y) - / dt G{zi;t)g{z2;t)G{z3, ..., Zk;y) 



y 
^G{zi;y)G{z2,...,Zk;y) - dt i—G{z2,zi;t)] G{z3, ...,Zk;y). 



Repeating this procedure one arrives at the following integration-by-parts identity: 

G(zi, ..., Zfe; y) + {-l)''G{zk, ..., Zi;y) 

= G{zi;y)G{z2, ..., zu; y) - G(z2, zi; 2/)G(z3, ..., Zk] y) + ... 

-{-lf-^G{zk-i, ...zi; 2/)G(zfc; y), (7.27) 

which relates the combination G(zi, ..., Zk\ y) + (— l)'^G(zfe, ..., zi; y) to G-functions 
of lower depth. This relation is useful in simplifying expressions. Eq. (j7.27|l can 
also be derived in a different way. In a Hopf algebra we have for any non-trivial 
element w the following relation involving the antipode: 

^u;(i) -5(^(2)) = 0. (7.28) 
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Here Sweedler's notation has been used. Sweedler's notation writes the coproduct 
of an element w as 

A{w) = ^u;(i)®u;(2). (7.29) 

(w) 

Working out the relation (|7.28|) for the shuffle algebra of the functions G(zi,..., 
Zkw), we recover (|7.27|l . 

We may now proceed and check if H7.28|l provides also a non-trivial relation for 
the quasi-shuffle algebra of Z-sums [99] . This requires first some notation: A com- 
position of a positive integer /c is a sequence I — (ii, ..., ii) of positive integers such 
that ii + ...ii = k. The set of all composition of k is denoted by C{k). Compositions 
act on Z-sums as 

{ii,...,ii) o Z{n-mi,...,mk;xi,...,Xk) 

= Z [n:mi + ... + 771^^,771^^+1 + ••• + "^n+is; •••7 "^-n+.-.+ii-i+i + •■• 

^"^^^l + ...^-■i^ ; Xi...Xi-^ , Xi-^^i. ..Xi-^j^i.^ ,...,X2j^ + ,..-|-2;_^-j-l..-^2 2^ + ...-j-i;j , yi .o\jj 

e.g. the first «i letters of the Z-sum are combined into one new letter, the next 
72 letters are combined into the second new letter, etc.. With this notation for 
compositions one obtains the following closed formula for the antipode in the quasi- 
shuffle algebra: 



SZ(n;mi,...,mk;xi,...,Xk) = (-1)' 2^ I o Z{n;mk, ■..,mi;Xk, ...,xi) 



\k 

iec(k} 



(7.31) 



From (|7.28(l we then obtain 

Z{n;mi, ...,mk\xi, ...,Xk) + {-if Z{n]mk, ...,mi;xk, ...,xi) 
= - ^ P^{Z{n\mi,...,mk]Xi,...,Xk)) 

adm. cuts 

■S {R'^{Z{n;mi, ..., 777fc; xi, ....Xk))) 

-(-!)'' Y^ IoZ{n;mk,...,mi;xk,...,Xi). (7.32) 

/ec(fe)\(i,i,...4) 

Again, the combination ^(77; 7771, ...,mk; xi, ...,Xk)+{—lfZ{n; irik, ..., mi;xk, --.jXi) 
reduces to Z-sums of lower depth, similar to (|7.27l) . We therefore obtained an 
"integration-by-parts" identity for objects, which don't have an integral represen- 
tation. We first observed, that for the G-functions, which have an integral rep- 
resentation, the integration-by-parts identities are equal to the identities obtained 
from the antipode. After this abstraction towards an algebraic formulation, one can 
translate these relations to cases, which only have the appropriate algebra structure, 
but not necessarily a concrete integral representation. As an example we have 

Z(77; 7771, 7772, 7773; xi,a;2,a;3) - Z(?7; 7773,7773, 7771; X3,X2,a;i) = 

Z(77; r77i; a;i)Z(r7; 7773, 7773; X2, X3) - Z(77; m2,mi;x2,Xi)Z{n; 7773; 0:3) 
—Z{n; mi + 7772; XiX2)Z{n; 7773; X3) + Z{n; 7772 -|- 7773, 7771; a;2X3, Xi) 
-|-Z(r7; 7773, nil + 1112; X3, 2:ia;2) + Z(77; r77i + 7772 + m3;xiX2X3), (7.33) 
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which expresses the combination of the two ^-sums of depth 3 as Z-sunis of lower 
depth. The analog example for the shuffle algebra of the G- function reads: 

G{zi, Z2, Z3;y) - G(z3, Z2, zi;y) = G{zi;y)G{z2, z^; y) - G{z2,zi;y)G{z^;y). 

(7.34) 

Multiple polylogarithms obey both the quasi-shuffle algebra and the shuffle algebra. 
Therefore we have for multiple polylogarithms two relations, which are in general 
independent. 

7.3 Numerical evaluation. At the end of the day of an analytic calculation 
physicists would like to get a number. This requires a method for the numerical 
evaluation of multiple polylogarithms. As an example I first discuss the numerical 
evaluation of the dilogarithm [49] : 

Li.(.) = -/rfti^^fj (7.35) 

"=i 

The power series expansion can be evaluated numerically, provided |a;| < 1. Using 
the functional equations 

Li2(x) = -Li2 Q) - ^ - i (ln(-x))% 

Li2(a;) = -Li2(f-a;) + — -ln(a;)ln(l-a;). (7.36) 

D 

any argument of the dilogarithm can be mapped into the region |a;| < 1 and — 1 < 
Re(a;) < 1/2. The numerical computation can be accelerated by using an expansion 
in [— ln(l — x)] and the Bernoulli numbers Bi: 

oo 

The generalisation to multiple polylogarithms proceeds along the same lines [100]: 
Using the integral representation 

Gmi,...,mfc {Zi, Z2, ••■, Zk\y) = (7.38) 

' 'dt \"^-' dt fdt \"^-' dt fdt V""--^ dt 



-o 



t J t — Zl \ t J t — Z2 \ t J t — Zk 



one transforms all arguments into a region, where one has a converging power series 
expansion: 



oo oo 



31 



Gmi,...,mk (Zl, ■■■, Zk;y) = /,■■■/, T- i \ ■ \mi 1 

1 ( y\^ 1 /y^^" 



(j2 + ...+,7fc)"'^ \Z2) {jkY"' \zk 



(7.39) 



The Art of Computing Loop Integrals 43 

The multiple polylogarithms satisfy the Holder convolution [66]. For zi 7^ 1 and 
Ztu 7^ this identity reads 

G(zi,...,z^;l)= (7.40) 

1\ „ / 1~ 



2_^{-iy G ( l-Zj, l-Zj_i,...,l-zi;l 1 G lzj+i,...,Zyj;- 

The Holder convolution can be used to accelerate the convergence for the series 
representation of the multiple polylogarithms. 

7.4 Related functions. In the literature one often encounters log-sine inte- 
grals. These are closely related to inverse binomial sums and can be expressed in 
terms of multiple polylogarithms. Here I briefly summarise the results from the 
literature. The following inverse binomial sum can be evaluated with elementary 
functions as follows: 

/n - r(n+2^x^ _ 2V^arcsin(Vi) 

In the literature, evaluations of inverse binomial sums of higher weights are given 
in terms of log-sine functions [101-108] : 

where 9 — 2 arcsin y/z and the log-sine functions are defined by 



Lsf\9) = - I d9' {9'f\n'-''-^ 



2 sin — 
2 



(7.43) 





By analytic continuation the log-sine functions arc then related to polylogarithms 
[108]. A simple example is given by 

Ls^°\9) - Ch{9), (7.44) 

involving the Clausen function CI2. The Clausen functions Cl„ are given in terms 
of polylogarithms by 

r ^ [Li„ (e*«) - Li„ (e-''«)] , n even, 
Cln{9) - <^ (7.45) 

I HLin(e'')+Li„(e-^»)] , n odd. 

8 Outlook 

In the previous section we discussed methods to perform the Feynman param- 
eter integrals. Unfortunately not every Feynman integral can be solved with these 
methods. In this section I will try to give an outlook towards open questions and 
future directions. One fundamental question is related to the type of functions 
which occur in the results of Feynman integrals. In sect. 18. il l briefly discuss the 
massless master two-loop two-point function. Here multiple zeta values are suffi- 
cient to express the result. However, the class of multiple polylogarithms (which 
contains the multiple zeta values) might not be sufficiently large to accommodate 
the results of all Feynman integrals. Indications come from the three-loop two- 
point function with equal internal masses, where complete elliptic integrals occur. 
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This is shortly discussed in sect. 18.21 In this article we considered mainly integrals 
associated to individual Feynnian graphs. The physical scattering amplitude is the 
sum over all Feynman graphs. It is not unusual that the result for a scattering 
amplitude can be expressed more elegantly than the results for individual Feyn- 
man graphs. In this context I review in sect. 18.31 a method for the computation of 
loop amplitudes based on unitarity, which by-passes traditional Feynman graphs. 
I conclude this section with a review of a conjecture, which states that the /-loop 
amplitude in maximally supersymmetric Yang-Mills theory is basically determined 
by the one-loop amplitude. 

8.1 The two-loop two-point function. In massless theories all two-loop 
two-point functions can be derived from the following master integral: 

f(2,5)/ \ -2/ 2\ 2^12345 -2m+2e /o i \ 

d^fci f d^k2 1 



i^D/2 J ,^D/2 (_)t2)-i (_^2)-2 (_^2)-3 (_^2)-4 (_^2)-5 ' 

where k^ = k2 — p, k/^ — ki ~- p, k^ — k2 — ki, D ~ 2m — 2e and 

_ r(i + £)r(i-£)^ 

'^ - — nr^) — ■ ^^-^^ 

The trivial dependence of this integral on (— p^) is already factored out in eq. H8.1|l . 
Therefore for given values of m and Vj , the Laurent series in e contains only num- 
bers. For a long time it has been an open question what type of numbers one 
encounters in this Laurent expansion. With the help of the methods discussed in 
sect. 16.41 it was possible to prove that only multiple zeta values occur in the Laurent 
expansion of the two-loop integral v'^'\ra~ 6,1^1, 1^2, 1^3, 1^4,1^5), if all powers of the 
propagators are of the form i>j — rij + ajS, where the rij are positive integers and 
the aj are non-negative real numbers [109]. As an example one has 

(1 - 2e) /(2>5) (2 - e, 1 + £, 1 + £, 1 + e, 1 + e, 1 + £) = 6C3 + 9(4^ + 372C5£^ (8.3) 
+ (915C6 - 864CI) e^ + (18450(7 - 2592(4(3) s'^ + (50259(8 - 76680(5(3 
-2592(6,2) £^ + (905368(9 - 200340(6(3 - 130572(5(4 + 66384(|) e^ + 0{e'^). 

8.2 Elliptic integrals. Laporta considered a three-loop two-point function 
with equal internal masses [110]. This integral is given by 

oc 00 

f dl f dm 

J A(/,-l,-l)i A(m,/,-l)A(m,-l,-l)- ^^'^^ 



Based on strong numerical evidence, this integral can be written as 

A^ K{w^)K{w+), w±^ ^^ z± = -(2-V3] U±VT5\ , (8.5) 

where K(x) is the complete elliptic integral of the first kind: 

1 
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8.3 The unitarity based method. We conclude the discussion of techniques 
for loop calculations with a method based on unitarity [6, 111]. I will discuss this 
method for one- loop amplitudes. Within the unitarity based approach one chooses 
first a basis of integral functions li E S. For one-loop calculations in massless QCD 
a possible set consists of scalar boxes, triangles and two-point functions. The loop 
amplitude A^^' is written as a linear combination of these functions 

.4(1) = E^^^»+^- (8.7) 

i 

The unknown coefficients c^ are to be determined, r is a rational function in the 
invariants, not proportional to any element of the basis of integral functions. The 
integral functions themselves are combinations of rational functions, logarithms, 
logarithms squared and dilogarithms. The latter three can develop imaginary parts 
in certain regions of phase space, for example 

Imln(^^^) = -n[e{s)~e{t)], 
ImLiall- [~'~'^} ] = -lnfl--)Imln 



{-t-iS)J \ tJ \-t-iS, 

Knowing the imaginary parts, one can reconstruct uniquely the corresponding inte- 
gral functions. In general there will be imaginary parts corresponding to different 
channels (e.g. to the different possibilities to cut a one-loop diagram into two 
parts). The imaginary part in one channel of a one- loop amplitude can be obtained 
via unitarity from a phase space integral over two tree-level amplitudes. With the 
help of the Cutkosky rules we have 

In^^''^ = ^^f rf^D,2^ -,,2^ -A ^Af. (8.9) 

J (27r)^ kf + iSk^ + tS ^ ^^ 

A'^^^ is the one- loop amplitude under consideration, A]^ and A^^' are tree- level 
amplitudes appearing on the left and right side of the cut in a given channel. 
Lifting eq. (|8.9I) one obtains 

^(1) ^ /(g7^^4"^4°^+ cut free pieces, (8.10) 

where "cut free pieces" denote contributions which do not develop an imaginary part 
in this particular channel. By evaluating the cut, one determines the coefficients Ci 
of the integral functions, which have an imaginary part in this channel. Iterating 
over all possible cuts, one finds all coefficients c^. One advantage of a cut-based 
calculation is that one starts with tree amplitudes on both sides of the cut, which 
are already sums of Feynman diagrams. Therefore cancellations and simplifications, 
which usually occur between various diagrams, can already be performed before we 
start the calculation of the loop amplitude. 

In general, a cut-based calculation leaves as ambiguity the ration piece r in eq. 
(|8.7() . which can not be obtained with this technique. One example for such an 
ambiguity would be 



{2tt)D k^{k~qY ' 



.11) 



This term does not have a cut and will therefore not be detected in a cut-based 
calculation. However, Bern, Dixon, Dunbar and Kosower [111] have proven the 
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following power counting criterion: If a one-loop amplitude has in some gauge a 
representation, in which all n-point loop integrals have at most n — 2 powers of the 
loop momentum in the numerator (with the exception of two-point integrals, which 
are allowed to have one power of the loop momentum in the numerator) , then the 
loop amplitude is uniquely determined by its cuts. This does not mean that the 
amplitude has no cut-free pieces, but rather that all cut-free pieces are associated 
with some integral functions. In particular A/" = 4 supersymmetric amplitudes 
satisfy the power-counting criterion above and are therefore cut-constructible. 

QCD does in general not satisfy the power-counting theorem and leaves as an 
ambiguity the rational function r. In principle one can obtain the rational piece r by 
calculating higher order terms in e within the cut-based method. At one-loop order 
an arbitrary scale p?'^ is introduced in order to keep the coupling dimensionless. 
In a massless theory the factor /i^^ is always accompanied by some kinematical 
invariant s~^ for dimensional reasons. If we write symbolically 

^'■' -~ hi3^' M"^' --i3^' ■ <-^) 

the cut-free pieces co(so//x^)^'^ can be detected at order e: 

co(^) ' - co-ecoln(^^)+0(e2). (8.13) 

8.4 Iterative structures of loop amplitudes. I would like to conclude this 
section with an outlook towards current research related to iterative structures of 
loop amplitudes. The loop amplitudes in maximally supersymmetric M — A Yang- 
Mills theory (MSYM) are a popular play-ground, as calculations in this theory 
tend to be simpler compared to the corresponding ones in QCD. In addition, one 
focusses on the leading-colour contributions to the loop amplitudes. Here only 
planar diagrams contribute. Let us denote by M„ = An /An the ratio of the 
l-loop amplitude to the corresponding tree amplitude. It is conjectured that the 
two-loop n-point function in MSYM is related to the one-loop n-point function 
by [112] 

M^^\e) = i(MW(e))' + /(2)(£)M«(2£)-^C4 + 0(e), 
/(2)(£) = -C2-C3e-C4e'. (8.14) 

This conjecture has been verified for n = 4. For n — b some partial results are 
available, which support the conjecture [113]. In addition it is believed that this 
iterative structure generahses to higher loops. For three loops it is conjectured [114] 
that 

M(3)(£) = -i(MW(e))' + MW(£)M(2)(£)+/(3)(,)MW(3£) 
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2i6+rO^^+rT+9 



7TT7 + n^i ) Ce + ( -- + 7TC2 ) Cl + 0{e) 



f^'He) = yC4 + (6C5 + 5C2C3)e+(ciC6 + C2C3)e' (8.15) 

Again, this conjecture has been verified for n = 4. The two constants Ci and C2 
cannot be determined from the n = A calculation. 
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9 Summary 

In this article I discussed loop integrals which occur in a perturbative ap- 
proach to quantum field theory. I reviewed standard techniques, which allow us 
to exchange the integration over the loop momenta against Feynman parameter 
integrals as wells as methods, which reduce tensor integrals to scalar integrals. The 
important sub-class of one-loop integrals was discussed in detail. The main part of 
this lecture was devoted to the computation of Feynman parameter integrals, with 
an emphasis on the mathematical structures underlying these computations. One 
encounters iterated structures as nested sums or iterated integrals, which form a 
Hopf algebra with a shuffle or quasi-shuffle product. In the final results multiple 
polylogarithms appear, and their properties have been discussed at length. The 
last section addressed open questions and conjectures. 
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